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Abstract 


This article deals with a special type of atmospherics produced within snow squalls. It 
discusses the possibilities of the genesis of such atmospherics. There is also a discussion of the 
theoretical principles and methods of measuring atmospherics by means of open-air antenna 
circuits and suitable cathode ray oscillographs. The article contains examples of typical meteoro- 
logical situations with atmospherics of the type discussed. Special attention has been given 
to a case where a bank of cumulo-nimbus clouds was located near the east coast of the 
Baltic and visible from Uppsala. Exemples are given of the typical variation forms of 

atmospherics produced in squalls within a limited range inside this cloud-bank. 


Introduction 


It is a generally accepted fact that terrestrial 
atmospherics are produced by lightning dis- 
charges from thunderstorm clouds. In such 
northern countries as Sweden this applies 
mainly to the warm season. Thunderstorms 
occur seldom during the cold seasons and they 
produce atmospherics of a variation feature 
typical of lightning discharges. On the other 
hand, we have observed atmospherics of such 
a number and intensity during the cold season 
in Sweden and when the sources of atmo- 
‚spherics within the range of measurement 
sensitivity could not have been influenced by 


thunderstorms that we are obliged to assume 
the existence of sources capable of producing 
atmospherics in other ways than by lightning 
discharges in the accepted sense of the word. 
This applies especially to those occasions when 
the general meteorological situation is charac- 
terized by the occurrence of snow squalls. 
Hence during the past winter we have under- 
taken some special investigations to ascertain 
whether or not electrical discharges within 
snow clouds, of other variation types than 
common lightning discharges, might serve as 
prominent sources of atmospherics. 
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Electric Field Conditions in Squalls of Snow 
Clouds. 


It is a well known fact that areas of dry-snow 
clouds, or ice crystals, may serve as sources of 
considerable static-charges. This is proved by 
experience when airplanes pass through arcas 
that produce heavy precipitation—static inter- 
ferences (Gunn, HALL, KINZER 1946 and GUNN 
1948). According to Gunn, the highest intensi- 
ties of the electrostatic field observed when 
aircraft pass through a thunderstorm-cloud 
were found at altitudes somewhat comparable 
to the freezing level. Mean field intensities of 
1,300 volts/cm were observed from nine 
different thunder-clouds. The free charge was 
characterized mainly by bipolar distributions; 
hence the field lines were distributed mainly 
within the cloud. 

Sometimes the intensively charged precipita- 
tion particles from a snow cloud squall must 
be distributed, e.g. in a volume extending 
from the cloud and down to the earth. If, as a 
first approximation, we assume this volume to 
be a cylinder and the volume charge distribution 
to be homogeneous, it is easy, by means of 
the potential theory, to calculate the field force 
intensity along the axis of the charged cylindric 
volume. This field intensity may be consider- 
able if the cylinder is high enough, as for in- 
stance in a high snow squall cloud. On the 
other hand, in ordinary thin snow clouds the 
electric field intensity will not be very great. 

I have personally observed such a case during 
a trip to U.S.A. in November 1928 with the 
Swedish liner Gripsholm. On a dark afternoon 
the ship was passing through a heavy snow 
squall, and on the masts and other pointed 
high parts of the ship I observed very pro- 
nounced Saint Elmo’s-fire phenomena. This 
was a typical evidence of the strong electrical 
field through which the liner was passing. 
There are many other indications of heavily 
charged particles in snow squalls. 

One evident and necessary condition for the 
production of electric discharges of a disrup- 
tive character, e.g. in the form of sparks, 
corona discharges, or streamers inside a snow 
cloud must be a field force of sufficient inten- 
sity. There is no evidence that the discharge 
resistance would be significantly diminished 
within an atmosphere filled with snow particles. 
On the other hand, the well known values be- 


tween spheres—30,000 volts/cm in ordinary 
laboratory air—might not be necessary, the 
more so if one considers the diminishing of the 
discharge resistance with decreasing atmos- 
pheric pressure. Values down to 3,000 volts/cm 
have, according to Gunn (1948), been measured 
directly just before a flying bomber was struck 
by lightning. The height of the bomber was 
about 4 km above ground. 

In this connection one particular comment 
appears appropriate. A cumulo-nimbus cloud 
as well as a snow squall cloud of about the same 
dimensions can be the site of considerable 
charges and hence of a strong electric field 
force between different regions of the cloud. 
If, in a cumulo-nimbus cloud, disruptive 
electric discharges in the form of lightning are 
observed, the cloud and the process will be 
classified as a thunderstorm from the meteoro- 
logical point of view. In other words, there 
exists a limiting stage of the intensity of the 
field force which is the necessary condition 
for the development of an electric discharge and 
hence for classifying the cloud as a thunder- 
storm. An ordinary flash of lightning between 
the cloud and the earth has spark lengths of 
between 1 and 2.5 kilometres; spark lengths 
between clouds reach much greater values. 
From an electrophysical point of view there 
must exist discharge phenomena of much 
shorter discharge lengths before such long 
discharges occur. In this connection we should 
keep in mind the existence of the typical 
predischarges connected with ordinary light- 
ning discharges. Such shorter discharges be- 
tween limited regions of the charged clouds 
must be assumed to have a low intensity of 
luminosity. They are, in general, hidden by 
the accumulation of particles. Thus they are 
not visible unless the observer is quite close to 
the discharge channel. An observer on the 
ground will not be able to notice these dis- 
charges. Our experience of the low visibility 
within snow squalls is sufficient to acknowledge 
this opinion. 

It is the author’s opinion that electric dis- 
charge of smaller range must occur within 
thunderstorm clouds and especially within 
snow squall clouds of a sufficient extension. 
Both spark lengths, current variation values, 
and hence transported charges are smaller 
than the values we have to accept in the case 
of ordinary lightning discharges, but the 
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theoretical considerations presented below 
show that atmospherics generated by the for- 
mer type of discharge do not differ in the 
general characteristics of their variation from 
what applies to a lightning discharge. The 
dimensions of field force and current variations 
of the case in question must evidently be much 
more limited. 

The only way to demonstrate the existence 
of the type of discharges mentioned is to use 
experimental electromagnetic methods. The 
intention is here to prove the existence of 
atmospherics caused by squalls of snow clouds. 
Here it is first necessary to present some 
theoretical considerations. 


Theoretical Considerations 


a) Effects of electric discharges in electric field 
above ground. 

It is quite easy to calculate the electromagne- 
tic variations resulting from an electric dis- 
charge of a cloud. To simplify the case we 
assume the discharge path towards the ground 
to be vertical and the transported charge to be 
concentrated to a point. As is well known, 
such an estimation applies to a homogeneous 
charge within a spherical volume—an approxi- 
mation which is acceptable in many cases. 

We assume a charged volume of the cloud 
to be concentrated to a point À at a distance d 
from o and at a height h above the ground 
(see Fig. 1). Further we assume that the trans- 
ported charge has values + Q, and the current 


value I at time £—— where c represents the 
c 


velocity of light. 

From the well known theory of electric 
images by Lord KELvIN we are able to consider 
the action of a discharge between point A 
and the ground as follows. Besides charge 


ERP 


Fig. r. Circuit for calculating variation of electromagnetic 
field from vertical lightning paths. 


+ Q in point A we have the same charge but 
with the opposite sign (— Q) in point A, 
symmetrically situated towards A with respect 
to the ground. In the discharge circuit two 
currents must be considered, both having the 
same intensity I, one passing from A to the 
ground and the other in the same direction 
from the ground towards Aj. 

From an analysis by LEJAY (1926) it is pos- 
sible to apply the general equations of Max- 
WELL to the discharge circuit. He has also 
given general solutions in the case of discharges 
towards the ground in other directions than 
the vertical. It is possible to carry out similar 
calculations also with respect to discharges 
between two regions inside the cloud. 

The electrical action in a point o on the 
earth, situated at distance d, can be represented 
by a vector F at right angles to the ground, 
and with a downward direction in the actual 
case. The expression for the vector in electro- 
magnetic units C.G.S. is given by equation (1) 
under the supposition that d is great in com- 
parison with h. 


2h; h 2h dI 
Fe bee pe 


Here, as before, c represents the velocity of 
light, or 3: 10!% cm/sec. By introducing 


— we obtain 


ir 


2h _ 2h dQ) ; 
BEEF Si ITB PAG 


If the first term in (1) is neglected, we ob- 
tain the magnetic field vector H, which is at 
right angles to the direction from o towards 


the path of lightning. The expression for H is 


_ 2h beh al (3) 
ica a : 


The three terms that enter the field force 
as given by (2) represent, in sequence, the 
electrostatic field, the induction field and the 
radiation field. In equation (2) the distance 

: jie er 
enters with a factor — in the first, with 7 


d? 


; : er 2 
in the second, and with- in the third term. 


d 


Hence the relative values of the terms vary 
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considerably with the distance from the chan- 
nel of lightning. This will be shown through 
some values of I, Q, and dI/dt which have been 
taken from some of our own researches (No- 
RINDER 1944). These values are given in table 1: 


There is also another reservation to be made. 
In equation (1) we have considered the dis- 
charge path to be vertical and to pass between 
the ground and a point A within the cloud. In 
the special case that we have to deal with here 


Table 1 


Electrostatic term 


Induction term 


Radiation term 


2 28e 2h I 2h dI 
c? — C= em 
Distance da a Ei Gh 
dI 
2 4 dt 
s coul 10 coul 10000 À 50000 A |500 A/usec |I 000 A/usec|2 000 A/usec 
5 km 14.4 V/cm 28.8 V/cm 0.48 V/cm 2.4 Vicm| :0.4 V/cm| os V/cm) 1.0. V/em 
Io » Tess 216) 002 rom > 0.6 » O2 00 CA OS 
food té 10 2» 13.6.0 >, 12 ro) Goa 0.02 » 0.04 » 0.08 » 
TOOOM DE ré Tome). (3:6. rom». in tom.» So dr) 0.002 } 0.004 » 0.008 » 


Table ı contains such values of I, Q, and 
dI/dt as we have been able to measure in actual 
cases of lightning discharges. 

When dealing with discharges within snow 
clouds we must assume much smaller magni- 
tudes of the corresponding values. In table 2 
the result is given of such a calculation. 


Table 2 
Electrostatic Induction Radiation 
term term term 
2 2h 2h I 2h dI 
istance Caters Oa mor 
2 a a d di 
à al 
| Q dt 
0.5 coul 1000 A 100 A/usec 
5 km 1.44 V/cm]4.8- 10 ”V/cm|8- ro ?V/cm 
10 » OS A) (Pew rote lafero 200 
200: TS TOM 1.200 Hy AR too” à 
soo anno u Imasıo >. las 


There are several reasons why one must be 
cautious in applying calculations based on (2). 
The main objection against a too far-reaching 
calculation is the rapid variations with time of 
the current that in most cases govern the elec- 
tric discharges of the atmosphere. This fact 
will cause a complex variation feature of dQ/dt 
and d?Q/dt? the result of which will be very 
complex forms of the calculated curves. This 
has been demonstrated by the present writer 
(NORINDER 1936). 


we must consider a discharge between two 
points within a cloud. Also in this case it is 
possible to use the theory of electric images, as 
has been done earlier by Wırson (1916) and 
applied by MAURAIN (1948). One finds from 
such calculations that either the higher or the 
lower charge will dominate the sign and 
magnitude of the electric force near the ground, 
depending upon the relation between heights 
of the charges above the ground and upon the 
distance to the point of observation. In the 
case of two equal charges of negative and 
positive signs respectively, one situated above 
the other, a change of the sign in the field 
force will occur at a certain distance from the 
discharge channel. If we assume a negative 
charge to be situated at a height of 2 km, and 
another, positive one, at a height of 5 km 
vertically above the other, the field force will 
change sign at a distance of 5.5 km. 

There are other considerations to be inlucded 
in the calculations. The discharge impulses 
pass through regions of the atmosphere where 
the ionization sometimes changes considerably 
from one locality to another. The discharges 
are characterized by very shifting variation 
features, which must be considered to be equi- 
valent to frequencies of different group velocity 
values. Hence the change of the ionization in 
the atmosphere will cause additional modifica- 
tions of the curve forms—results that are not 
possible to calculate. A further complication 
arises from the effect of the reflection of the 
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impulses at the upper layers of the ionosphere. 

The calculations presented above will be 
considered as a help in understanding the trans- 
formation that an electric discharge impulse 
undergoes with regard to distance. At the 
source within the cloud the impulse current 
variation may have a typical and unipolar 
exponential variation form with rapid super- 
imposed shorter variations. When the impulse 
has passed a distance of some kilometres it has 
very often undergone a transformation into 
variation type of double polarity. During their 
subsequent passage the transformations will be 
much more complicated. 

The way to carry out research into the 
interesting transformation of the impulses with 
regard to distance must primarily be a problem 
reserved for systematic experimental measure- 
ments. For research of this kind it is necessary 
to analyze the transformations of the atmos- 
pherics in individual cases by using recording 
stations at known distances from the sources 
of the discharges. One condition is that measur- 
ing methods are available which allow results 
reducible to an equivalent measure. 


b) Methods of measuring electric field variations 
at distances from discharges. 


The method of measuring the variations of 
the electric field components as emitted by 
a discharge at various distances from a current 
path is based on the application of a suitable 
open antenna circuit. This arrangement was 
used by the present writer (NORINDER 1924, 
1925, 1934) when the aperiodic variation forms 
of lightning discharges were first discovered. 
The same method was applied by APPLETON 
and WATT (1923) in the first study of the 
variation forms of atmospherics. 

The antenna circuit must be placed parallel 
with the horizontal level of the ground in an 
open and undisturbed region (see Fig. 2). Hence 
we are able to consider the antenna wire to be 
exposed mainly in an equipotential surface. 
Let us consider the field force to be E in volts/ 
metre; thus the potential above the ground 
at the height h is hE. In this case we suppose 
the field force lines to be directed down to the 
ground. In order to eliminate the oscillations 
proper of the antenna circuit provoked by 
some rapid variations of the lightning dis- 
charge, a damping resistance R, is introduced 
as well as a measuring resistance R,. Above 
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this latter resistance voltage variation V is 
connected to the deviation plates of an oscillo- 
graph. Parallel with the deviation plates a 
variable capacity can, if necessary, be intro- 
duced in order to regulate the voltage devia- 
tion. If this capacity is introduced it must be 
considered in the following calculations. 

The relation between V and E is given by 
equation (4) 


dE R, + R,fdV I A, 
Te nee a 


We have in a special case R, < R, (e. g) 
R, = 2502, R, = 20000), and :R, can thus 
be omitted in comparison with R,. Taking this 
into consideration, and with R, = R, we obtain: 


an a 
HI ae Re (s) 


The sign of the right hand member depends 
on how the sign of the oscillographic voltage 
is determined. By using the sign of E, so that 
the lines of force go towards the ground and 
by supposing an increase of the field force 
with time, the result will be a positive sign of 
the right member of (5). Which of the two 
terms dV/dt or V/RC is in preponderance 
depends on the sloping ratio of the recorded 
curve and on the value of RC. With a small 
value of R, we have approximately: 


dE I I 


In this case the second term is the prepon- 
derant one. In order to obtain the variation of 


E 


Fig. 2. Schematic diagram of antenna circuit for 
calculation of field variations of lightning discharges 
and atmospherics. 
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field force E, we must in such a case integrate 
the oscillographic curves. This is called the 
dE/dt method. With a high value of R, (s) 
will be transformed to: 


dE tadV I 
as Bay ak es een 


In this case the first term is the more im- 

ortant one and we receive the variation of 
the field force without an integration proce- 
dure and thus directly. This is called the E 
method. 

The more elaborate dE/dt method is used in 
those cases when the field measurements are 
taken at a place where local disturbances, e. g. 
from high tension line systems, disturb the 
results. It is always preferable to select a locality 
for field measurements where no local disturb- 
ances can influence the measurements. In 
such a case one has ample opportunity of 
using, if necessary, both measuring methods. 
With the exception of measurements made 
quite near to the discharge channel it is neces- 
sary to introduce into the measuring circuit 
an amplifier of such a construction that the 
atmospherics can be recorded without dis- 
tortion. Such a circuit with an amplifier in- 
cluded is reproduced in Fig. 3. 


Fig. 3. Antenna circuit connected above amplifier to 

cathode ray oscillograph. Constants used: Ca = 300 uuF: 

Cc = 0.001 UF; Cp = 0.001 — 0.1 MF; Rp = 5 meg Q; 
Rs = 50 (2. 


It is necessary to operate the field station at 
places where local disturbances are eliminated. 
An important condition must be, that the field 
of the antenna circuits are undisturbed and 
commensurable. Such a field station, Funbo- 
Lövsta, one of the stations belonging to the 
Institute, is illustrated in Fig. 4. 

As a suitable instrument for the measure- 
ments a special type of cathode ray oscillograph 
was developed. The cathode ray tube was of 
the sealed-off pattern. The cathode ray oscillo- 


Fig. 4. External view of a field station at Funbo-Lövsta 
for measuring lightning discharges and atmospherics. 


graphs constructed at the Institute were spe- 
cially adapted for use in the study of atmos- 
pherics. The instruments were provided with 
a film camera arrangement in combination 
with a sensitive lens system. This resulted in 
very good photographing properties also in 
the case of rapid phenomena down to a dura- 
tion of a microsecond. A view of such an 
instrument is reproduced in Fig. $—6. 


Fig. 5. External view of a 

cathode ray oscillograph 

with a tube of the sealed- 
off pattern, 


Fig. 6. Internal view of cat- 
hode ray oscillograph in 
Figs ss 
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The sealed-off type 
of cathode ray oscillo- 
graph has been used 
especially in simultan- 
cous recordings of at- 
mospherics at two sta- 
tions located at con- 
siderable distances from 
each other. In such cases 
the oscillographs were 
provided with a simul- 
taneous time marking 
system, controlled by 
syncronisation from the 
Swedish  syncronised 
electric power system. 
An original photograph 
of oscillograms of at- 
mospherics is repro- 
duced in Fig. 7. 

In order to localize 
the sources ofthe atmos- 
pherics it was necessary 
to operate with direc- 
tion finders. A view of 
such a direction finder 
constructed by STOFF- 
REGEN (1947) and de- 
scribed elsewhere is 
given in Fig. 8. 


Measurements of 
Atmospherics with 
Special Reference to 
Snow Squalls 


The most significant 
variation feature of the 
atmospherics emitted 
during the winter from 
snow squalls is their 


strikingly regular type. 


A 


Fig. 7. Original photo- 
graph of oscillograms of 


atmospherics. à : : 
observation this varia- 


tion type may be repeated for many hundreds 
of times within a very short time interval. 
By comparison one finds that the correspond- 
ing variation features of atmospherics emitted 
from thunderstorms at corresponding distances 
are characterized by a much more pronounced 
irregularity. That this must be the case follows 
from the shifting variation types which always 


During one period of 


Fig. 8. Direction finder specially constructed for locating 
sources of atmospherics. 


characterize lightning discharges. Also in cases 
where the sources of lightning discharge have 
been situated at considerably greater distances 
than those within the squalls of snow clouds 
we find a more pronounced irregularity in 
those atmospherics which are a consequence 
of flashes of lightning. 

It must be remarked that in all the following 
measurements the E method has been used for 
the vertical component of the electric field. 

In order to illustrate the difference referred 
to above, a few oscillograms, obtained from 
the two stations Funbo-Lövsta (Uppsala) and 
Orsa, at a distance of about 300 km from each 
other, are reproduced in Fig. 9—10. Simulta- 
neous records of atmospherics were taken at 
these two stations. Fig. 9 represents atmos- 
pherics of the long duration type. In this case 
the source was situated somewhere in Central 
Europe at a distance from the Uppsala station 
of about 1,500 km. Similar variation types from 
snow squalls have never been observed. In Fig. 
ro we have instances of another case with a 
distance of about 500 km from the source 
to the Uppsala station. The source was located 
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Fig. 9. Simultaneous records of atmospherics at two sta- 
tions separated by a distance of 300 km. 
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Fig. 10. Simultaneous records of atmospherics at two 
stations separated by a distance of 300 km. 


in the southern part of Sweden. If we compare 
these oscillograms of typical irregularities with 
the following oscillograms from snow squalls 
having regular variation forms, we shall find 
a very prominent difference. 

We shall now illustrate a situation which is 
very interesting with regard to the numerous 
repetitions of atmospherics of approximately 


similar type. This occurred on November ıst, 
1948—a day characterized by a cold front 
located over the North Sea. During this day 
we observed a very intense production of 
atmospherics from the front. Most of them 
were of such a regularity and variation form 
as to exclude the possibility that they might 
be caused by lightning discharges. Of about 
300 atmospherics recorded during a continuous 
period of only ten minutes, the ones illustrated 
by Fig. 11a, b show variation types of such 
a resemblance as to give the impression that 
all have been produced by a similar process in 
a special regular way. 

We shall give an instance of another case. 
On the ıst of March 1949 a snow storm was 
sweeping over Denmark, the northern part 
of Germany and on into Poland. Recording 
with a direction finder we observed a great 
number of atmospherics emanating from the 
directions mentioned. In Fig. 12 a, b some 
layed-out oscillograms are reproduced. There 
is not the least impression that any of them 
could have been generated by lightning dis- 
charges. They are typical of snow squalls. 
The measured values of the field forces agree 
well in their order of magnitude with the 
calculated field force values of table 2 at dis- 
tances of 1,000 kilometres. 

During the past winter season we have been 
on the look-out for a situation in which a more 
restricted region of snow squalls appeared 
not too far away from the observation station. 
As a further condition it was stipulated that 
no other sources of disturbances should inter- 
fere by producing atmospherics of the same 
or a higher order of magnitude. 

Such a situation occurred on the 22nd of 
November, 1948. On this day and the next we 
observed from the Uppsala station a stable and 
characteristic bank of cumulo-nimbus clouds 
clearly visible between 50° and 190° from the 
field station Funbo-Lévsta at a distance of 
10 kilometres from Uppsala. A picture of this 
cloud-bank is reproduced in Fig. 13. The situa- 
tion was observed to be very appropriate for 
carrying out investigations with the aim of 
ascertaining whether the stable cloud-bank of 
a cumulo-nimbus type was producing atmos- 
pherics of the special type mentioned above. 

Hence, at Funbo-Lövsta measurements of 
the variation of the electric field by the above 
mentioned E method were started. Efforts 
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Fig. 11 a, b. Characteristic types of 300 atmospherics of similar variation type which emanated from fronts in the 
North-Sea during a period of 10 minutes. 


were made to concentrate the records on such 
types of atmospherics as emanated from a 
certain region of the cloud-bank. This was 
possible by simultaneous measurements of 
the source of the atmospherics with the direc- 
tion finder. j 
During the observation periods of the day 
in question the sources of the atmospherics of 
the special variation form and characterized 
by a high intensity level were located at direc- 
tions from 160° to 175° and, in some few cases, 
up to 180°. The level of the field force varia- 
tion of the atmospherics was very suitable in 
excluding records of atmospherics from other 
sources and distances. The latter could thus 
be pressed down to the zero-band of the 
oscillograph. The amplification was regulated 
in such a way that about 20 % of the recorded 
atmospherics resulted in a full deviation on the 


fluorescent screen of the oscillograph. It was 
not possible to avoid that some of the atmos- 
pherics became somewhat oversteered. This 
was necessary in order to permit recording 
also of the small amplitudes of the predis- 
charges. In the case of lower amplitudes the 
predischarges were covered within the zero- 
band of the oscillograph, and thus the record 
begins with the main discharge. To avoid this 
difficulty it would be necessary to operate 
simultaneously two oscillographs regulated at 
different levels of sensitivity. Reproductions 
of the original records are given in Fig. 14—15. 
In Fig. 14 the amplitude level was too high, 
which was necessary to obtain the predis- 
charges in question. 

Some examples of layed-out records of the 
atmospherics are reproduced in Fig. 16 where 
two atmospherics are shown having a typical 
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Fig. 12 a, b. Atmospherics from snow-storm front. 


Fig. 13. À bank of cumulo-nimbus clouds over the Baltic Sea, photographed from the Institute. 
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Fig. 14. Original oscillograms of atmospherics from 
snow-storm cloud. 
Original oscillograms of atmospherics from 
snow-storm cloud. 

Fig. 16. Layed-out records of atmospherics from snow- 
storm cloud with predischarges. 

Layed-out record of atmospherics from snow- 
storm cloud of long duration. 
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initiating period with small amplitudes and 
lasting about 250 usec. This period is followed 
by the main discharge with high amplitudes 
and regular variation forms. The total dura- 
tion attains about 1,000 usec. An interesting 
variation type of the rather long duration of 
approximately 1,750 psec. is reproduced in 
Fig. 17. In this case we have two discharges in 
the same channel. Examples of the charac- 
teristic regular variation form of the atmos- 
pherics are shown in Fig. 18—19. If we compare 
this regular variation form with the charac- 
teristic lightning discharge variations of Fig. 
9—10, the marked difference becomes evident. 

In order to obtain a general survey of the 
weather situation during the day of observa- 
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Fig. 18. Atmospherics of regular variation form from 
snow-storm cloud. 


tion we have from the Swedish Meteorological 
and Hydrological Institute received the obser- 
vations from some of the nearly meteorological 
stations on the Swedish east coast and from 
adjacent stations on the coast of Finland. It is 
to be noted that none of the stations reported 
lightning discharges. Some of them, especially 
Ut, reported snow squalls at the observation 
periods 07%oo and r9oo, and during interven- 
ing observation periods this station reported 
cumulo-nimbus clouds and a cloudiness of 
5/19 —"I1o- 
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Fig. 19. Atmospherics of regular variation form from 
snow-storm cloud. 


THE MEASUREMENT AND NATURE OF ATMOSPHERICS 13 


From an examination of the values received 
by the direction finder it follows that the source 
of the recorded atmospherics was situated 
within a comparatively narrow region not 
too far away from the Swedish meteorological 
station Utö. The distance from the Uppsala 
station is estimated to be at the order of 120 to 
150 kilometres. These distances agree rather 
well with the amplitude values of table 2, 
from which follow amplitudes of the ver- 
tical field force dE/dh of the order of 0.5 
volts/metre. 

The following general conclusions appear 
to be acceptable: As a sometimes important 
source of atmospherics during winter seasons 
in countries of higher latitudes clouds of the 
cumulo-nimbus type bringing heavy snow 
squalls must be taken into consideration. From 
different points of view, it seems evident 
that such squalls must be the site of considerable 
electric charges per unit of volume and extend- 
ing over considerable volumes. There is no 


reason why the intense electric fields arising 
from these charges should not produce disrup- 
tive electric discharges. Compared with the 
ordinary lightning discharge channels the dis- 
charges within snow squalls are of a lower 
intensity and a shorter spark length. But there 
is no physical reason why these discharges 
should not, in spite of their weaker character, 
be intense enough to produce typical electric 
disturbances in the form of atmospherics. 


Evidently the first experiments dealt with 
here, in which a method of localizing the emis- 
sion of atmospherics from snow squalls was 
used, must be completed by simultaneous 
records from stations located at varying dis- 
tances from the source of emission. It must be 
of a special importance to try to operate in the 
vicinity of the source. It is also necessary to 
localize the site of separate discharges by a 
simultaneous operation using two direction 
finders along a suitable base line. 
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Abstract 


This paper contains an aerological analysis of a series of rapid-moving frontal waves 
associated with a well developed zonal current over central Europe. The gradual de- 


struction of this zonal flow through a retrograde blocking “wave’’ is 


described. The 


relationship between the blocking process and the deepening of series of waves ap- 
proaching the blocking zone from North America is discussed. The observational data 
and conclusions from this aerological study are presented for the purpose of rendering some 
assistance to theoreticians investigating atmospheric wave motions, by providing a few nu- 
merical values for certain characteristic parameters of these waves and by calling attention to 
a factor which appears to play a significant röle in the deepening of certain types of 
extra-tropical wave cyclones, viz. the variation with longitude of the basic current pattern. 


Introduction. 


The objectives of the researches described 
below are threefold. The well-known investiga- 
tions conducted by WILLETT (1944), NAMIAS 
(1947) and others during the last ten or twelve 
years have shown that during the winter 
season the general atmospheric circulation in 


middle de has a tendency to fluctuate 
between two extreme states, characterized by 
the predominance of either zonal or meridional 
motion. The successive stages in the break- 
down of a well-defined zonal circulation into 
a complicated pattern of meridional currents 
has been described on several occasions, but 
mostly with the aid of consecutive overlapping 
mean charts in which time averages of the sea 
level pressure distribution during five or seven 
day intervals have been plotted. In the present 
paper a particularly striking case of such a 
break-down of a fully developed west wind 
belt into a series of quasi-stationary cyclonic 
and anticyclonic vortices will be described 
with the aid of ordinary synoptic charts for 
the soo mb level and by means of other 


synoptic tools. As a by-product the case in- 
vestigated provides some information con- 
cerning the manifestations and westward 
speed of a typical ‘‘blocking wave”. Of late, 
the great importance of this process has been 
recognized by synoptic and theoretical me- 
teorologists alike but as yet no adequate 
synoptic description has been published. 

The second objective has been to provide a 
few empirical data for determination of the 
numerical values of the characteristic parame- 
ters of typical fast-moving frontal waves. 
The total number of investigations of polar 
front disturbances published since J. BJERKNES’ 
epoch-making paper appeared in 1919 is 
certainly very great. However, in the light of 
recent experiences it appears safe to conclude 
that this collection of case histories is extremely 
heterogeneous, including typical frontal waves 
as well as dynamically quite dissimilar major 
storms associated with the deepening of pla- 
netary wave-troughs in the upper west-wind 
belt. Furthermore, in most cases the systems 
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investigated have been treated as solitary 
phenomena; with such a starting point one 
is not likely to arrive at numerical data suitable 
to serve as an adequate basis for a realistic 
wave theory. The present study of wave cyc- 
lones makes no claim to pry as deeply into 
the structure of the individual wave cyclone 
and its weather manifestations as do some of 
the erlier studies by BERGERON, Swosopa 
(1924) and others, but a deliberate effort has 
been made to bring out those factors which 
from our present point of view appear likely 
to be of importance in the dynamic inter- 
pretation of these waves. 

The third objective of our studies has been 
to shed some light on the nature of the frontal 
waves and particularly on the occlusion 
process. Until recently, it has been assumed 
that the important changes in the flow pattern 
of the upper troposphere, in particular the 
break-down of the zonal circulation, are pro- 
duced by the occlusion of frontal waves. 
While it is true that a marked deepening and 
cutting-off of an upper trough in most, if not 
all cases, is accompanied by a marked deforma- 
tion of the fronts at the ground in a manner 
which perhaps best is described as an occlusion 
process (see Fig. 34 in NYBERG 1945), cases 
have been observed in which a whole series 
of disturbances of the short and rapid-moving 
wave type occlude without bringing about 
significant modifications of the upper flow 
pattern. Such a series of occluding waves is 
described below. It will be shown that the 
gradual retardation of the occluding waves is 
associated with a pronounced eastward 
decrease in the strength of the zonal ‘‘steering”’ 
current. One is thus led to two conclusions, 
one definite and one tentative: first, that the 
process of occlusion by itself is insufficient to 
bring about significant changes in the upper 
flow pattern and secondly, that the occlusion 
of fast-moving frontal waves just possibly may 
be the consequence of the rapid decrease 
downstream of the basic current in which 
they are embedded. If this view point is correct 
the occlusion process would be comparable 
to the amplitude increase and breaking of 
ordinary sea surface waves which occurs 
whenever the phase velocity of these waves is 
reduced, by decreasing water depths, below a 
prescribed critical value. The occlusion would 
thus represent an instability of form only; 


the total energy available to bring about 
these form changes would merely be that of 
the wave motion itself. This situation differs 
fundamentally from the one prevailing in 
those unstable waves which are capable of 
withdrawing potential and kinetic energy 
from the basic barocline current system in 
which the waves are embedded. 

In a further attempt to shed some light on 
the nature of frontal waves a comparison has 
been made between the family of frontal 
waves described below and a group of fast- 
moving waves travelling eastward across the 
Rocky Mountains in the northwestern part of 
the United States and studied by Hess and 
WAGNER (1948). The comparison brings out 
a striking agreement in the upper-level char- 
acteristics of these two series of disturbances, 
though the surface manifestations of the waves 
for obvious reasons are quite dissimilar. 
Further investigations may be required before 
final and definite conclusions can be drawn 
but the comparison nevertheless suggests that 
the short progressive frontal waves observed 
in the lower tropospheric layers, in spite of 
their obvious primary importance in as far as 
weather phenomena go, must be regarded as 
induced manifestations of relatively short, 
rapid-moving upper level waves associated 
in an as yet unknown manner with the jet 
stream zone and its strongly developed so- 


lenoidal field. 


I. Zonal Circulation and Frontal Waves. 
By R. BERGGREN. 


During the second week of February 1948 
the weather in western and central Europe 
and over southern Scandinavia was controlled 
by a series of fast-moving frontal waves. 
Large amounts of precipitation and high mean 
temperatures during the week were observed 
at most of the stations. The waves moved 
from the Atlantic over the British Isles and 
down to the Black Sea with a great regularity. 
The pressure centre and the occlusion point 
for each disturbance kept within a belt of 
about $00—600 km width. The narrowness of 
this strip (shadowed in Fig. 1) indicates the 
similarity between the paths of the individual 
disturbances. Figure 2 demonstrates the pe- 


Fig. 1. 


Cyclone track belt for February 6—12, 1948. 


riodicity of the waves. The slanting lines, 
which represent the consecutive positions 
of the wave centers, have approximately the 
same inclination, i.e. the translational velocity 
from west to east was nearly constant from 
wave to wave, which could not have been 
the case unless the upper, steering current 
remained nearly stationary throughout the 
period. The positions of the wave centers in 
this figure were determined from the 0700 
GMT positions of the occlusion points. The 
diagram as reproduced includes five of a total 
of seven well defined waves observed before 
the basic current pattern broke down. With 
the aid of this diagram it is possible to compute 
an average period for the frontal waves. One 
finds for this period 31 hours at 30° W, 32 
hours at 0°, and 35 hours at 20° E, and for the 
entire region an average period of 33 hours. 
The average wave speed over the Atlantic 
and Western Europe was about 70 km p.h. 
and the average wave length in the same 
region 2,500 km. Over Central and Eastern 
Europe the waves were retarded, as indicated 
by the break in the lines of constant phase. 

In this region the wave length was reduced 
to about 1 600 km and the speed of propaga- 


1 The westward motion of the break in the curves 
of constant phase in Fig. 2 may be taken as an indication 
of a westward displacement (retrogression) of a decrease 
in the zonal circulation. This retrogression may possibly 
be related to the “blocking action” described by Namıas 
and Crapp (1944). The speed of the westward motion, 
7—8 degrees of longitude per day at 50° N, is in good 
agreement with the results obtained by these authors. 
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tion to so km p.h. If we compute the fre- 
quency, or the ratio c/L, where c is the speed 
of propagation and L the wave length, one 
obtains, for the western part of the tracks, 
2.8: 10-2 hours 1, and for the eastern part 
3.2: 10°? hours-!. As the values for c and L 
are smoothed values it appears permissible 
to conclude that the frequencies (and periods) 
are approximately constant. 

To demonstrate the quasi-stationary character 
of the upper ‘‘steering” current during the 
period under investigation a mean chart 
(Fig. 3) has been computed from the seven 
individual soo mb morning charts for the 
time interval February 6—12, inclusive. This 
chart may be considered as typical of those 
cases in which the weather over western 
Europe and Scandinavia is controlled by 
frontal waves coming in from the west or 
southwest. The (sea level) frontal system in 
Fig. 3 is taken from the sea level chart for the 
morning of February 8. The general orienta- 
tion of this frontal zone agrees quite well with 
that of the mean current at 500 mb. Also 
this chart shows that the distance between 
consecutive waves decreases eastward. If the 
mean chart in Fig. 3 is compared with the 
individual charts for the same period one 
finds that the changes occurring during this 
week are of a minor nature. As an example 
of an individual chart we have reproduced the 
morning chatt for February 9, in the middle 
of the period (Fig. 4 a). There is good agree- 
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Fig. 2. Phase diagram for February 6—12, 1948, showing 
successive positions of the occlusion points of five con- 
secutive frontal waves. 
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Fig. 3. Mean soo mb chart for February 6—12, 1948. The frontal system is taken from the sea level chart for 
February 8, 0600 GMT. Contour heights given in geodynamic decameters. 


ment as far as the position of the main ridges 
and troughs are concerned. Also with respect 
to the steering current there are only minor 
differences between the two figures. 

Fig. 4 a brings out certain other noteworthy 
features. A comparison with the sea level 
chart for the same day (Fig. 4 b) shows that 
the fairly intense wave cyclone observed 
that day over the North Sea was reduced to a 
minor trough over Scotland at the soo mb 
level (it should be noted that there is a slight 
difference in time between the two charts). 

Fig. 4a shows that at the soo mb level 
extremely cold air had accumulated in the 
upper level trough east of Hudson Bay. On 
this, as well as on several subsequent days, an 


isotherm for — 50° C can be traced with the 
aid of the radiosonde data from the north- 
eastern part of North America. It will be 
shown in the next section that this cold air 
plays an important röle in the subsequent 
break-down of the zonal circulation. 

The sea level chart in Fig. 4 b is typical of 
the entire period. The well developed high 
over Russia with a ridge over central and 
southern Scandinavia is important for later 
events. It is this system which, a weck later, 
brings down very cold air over Central 
Europe. 

An examination of the variations of the 
geopotential of the main isobaric surfaces 
shows that the wave disturbances which 
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Fig. 4a. soo mb chart for February 9, 0300 GMT. Heights in geodynamic decameters. 


passed Great Britain during this period can 
not be considered as shallow phenomena. 
Fig. 5 gives the geopotential variations of the 
300, 500 and 700 mb surfaces (and the varia- 
tions of the sea level pressure) over the radio- 
sonde station Aldergrove in northern Ireland. 
The curves demonstrate that the individual 
waves extended through the entire tropo- 
sphere, though the particular disturbance that 
passed Aldergrove on February 7, 0900 GMT 
was considerably weaker than the others and 
did not acquire full intensity until it reached 
Germany. In spite of the great depth of these | 
disturbances and in spite of the fact that all of‘ 
them occluded they did not bring about any signift- | 
cant major changes in the upper flow pattern. | 
Fig. 4 b. Sea level chart for February 9, 1948, 0600 GMT. Aerological experience indicates that this re- | 
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Fig. 5. Diagram showing the height variations of the 

main isobaric surfaces, 300 mb, soo mb, and 700 mb. 

P, is the mean sea-level pressure in mb. Ordinate inter- 

vals in units of 200 feet, and for the sea level pressure 
10 mb. 
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Fig. 6a. Cross section for Feb. 12, 1500 GMT, 1948, 
showing isentropes (dotted lines; potential temperature 
in °K) and isopleths for observed wind speed (full lines; 
speed in m.p.s.) perpendicular to the section. The heavy 
line represent the frontal boundaries and the tropopause. 
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sult is characteristic of fast-moving frontal 
waves approaching Europe from the West. 
A series of N—S cross sections have been 
drawn, mainly through the British Isles. Fig. 
6a and 6b contain a section for February 12 
from Cornwall to the Shetland Islands with a 
characteristic picture of a strong westerly air 
current. Fig. 6 a shows isentropes (dotted lines) 
and isopleths (full lines) for the observed wind 
speed perpendicular to the section, i.e. re- 
presenting a wind direction of 290°. Fig. 6b 
shows isentropes (full lines) and isotherms 
(dotted lines). The wind maximum (jet 
stream) is situated at the 250 mb level, im- 
mediately beneath the tropopause, with the 
strongest vertical increase of the wind observed 
between the frontal layer and the 400 mb 
level. In this particular case the frontal layer is 
well-defined and seems to extend into the 
lower stratosphere over the Shetlands. Over 
southern England there is a weak subsidence 
inversion. A strong horizontal north-south 
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Fig. 6b. Cross section for Feb. 12, 1948, 1500 GMT, 

showing isentropes (full lines; potential temperature in 

°K) and isotherms (dotted lines; temperature in °C). The 

heavy lines represent the frontal boundaries and the 
tropopause. 
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Fig. 7. Mean cross section for Feb. 6—12, 1948, showing 
isentropes (dotted lines; potential temperature in °K) 
and isopleths for computed geostrophic wind speed 
(full lines; speed in m.p.s.) perpendicular to the section. 


temperature gradient exists between 52° N 
and 58° N. The figure shows that the solenoids 
are not concentrated exclusively to the frontal 
layer but that they are strongly developed 
also in the warm air. 

The mean cross section from Thorshavn on 
the Faroes to Cornwall (Fig. 7) gives an in- 
teresting picture of the mean distribution 
during this period of the westerly wind 
(computed with the aid of the geostrophic 
assumption.) A rather well-developed wind 
maximum is formed a little farther to the 
north than in the individual section in Fig. 6. 
This mean maximum is situated on the central 
line of the trajectory belt in Fig. 1. 

During the period under investigation the 
jet stream undergoes an interesting variation 
which appears as a fluctuation in the direction 
and strength of the current as well as in the 
height at which its center occurs. When a 
trough passes a particular station the wind 
in the jet shifts from southwest to northwest. 
At that time its speed reaches values from 50 %, 
to 100 % greater than between troughs. The 
tropopause sinks to its lowest height as the 
trough passes and rises rapidly afterwards. 
At the same time the height of the jet stream 
center increases. This is indicated by the slant- 
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ing isopleths at higher levels in Fig. 8 and by 
Table? 


Table 1. Variation of the height of the jet stream 
center with time, as indicated by the observed 
wind speeds (in knots) over Larkhill (51° N, 2° W). 


February 10 | February 12 


Time (GMT)| 09 dame ||P al 09 | 15 21 
200 mb| — 69 | 51 — 70 83 
250 70 548193 — | 110 | 109 
300 82 105 | 90 | 100 | 112 96 
350 87 114 69 109 QI 89 
400 89 104 | 60 | 102 OI 76 


The variation in the height of the jet stream 
center stands out quite clearly in this table. 
The fact that the table makes use of winds 
at constant pressure levels does not alter the 
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Fig. 8. The variations with time of observed wind speed 
(in knots) and wind direction over Larkhill (51° N, 2° W). 
M stands for maximum and m for minimum wind 
speed. The meaning of the different shadings in the wind 
direction diagram is given in the wind rose. 
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Fig. 9. Diagram showing the zonal wind speed (computed 
geostrophically) for different isobaric surfaces over an 
area of the same order of magnitude as a specially selected 
typical frontal wave. The area moves with the wave so 
as to occupy the same position relative to the disturbance 
at all times. The broken line represents the eastward 
velocity of the wave as determined from the twelve 
hourly displacements, (On February 10 the material for 
the 300 mb analysis was sparse, but the indicated drop 
of the geostrophic wind speed is beyond all doubt real.) 


result significantly since the greatest shift of 
an isobaric surface at Larkhill during these two 
days is of the order of magnitude 60 m in 6 
hours while the shift of the jet stream center 
is of the order of magnitude 1 km per 6 hours. 
The secondary maximum on February 10, 
1500 GMT at 200 mb is connected with the 
upper part of a double tropopause. An examina- 
tion of the jet stream thus gives the result 
that it undergoes great variations as the fast- 
moving frontal waves pass by both with 
regard to strength, direction and height. 

In order to investigate the connection be- 
tween the velocity of the disturbances and 
the upper current the diagram reproduced in 
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Fig. 9 was constructed. A moving, quasi- 
rectangular area covering ten degrees of la- 
titude and thirty degrees of longitude was 
chosen every twelve hours, in such a manner 
that the particular disturbance which on Feb. 9, 
0700 GMT (Fig. 4b) was situated over the 
North Sea always remained in the center of 
the area. The zonal wind component for this 
strip was then computed geostrophically and 
plotted in the diagram which gives the zonal 
wind components for the 1,000, 800, 700, 500 
and 300 mb surfaces during the period Feb. 
8, 0300 GMT to Feb. 10, 0300 GMT inclu- 
sive. The broken line gives the eastward 
velocity of the disturbance measured as a 
twelve hour mean and defined as the mean 
value of the displacement in longitude of 
the pressure center and of the occlusion point. 
The diagram brings out, in a striking fashion, 
two facts: firstly, the disturbance moves faster 
than the zonal wind up to a certain critical 
height at which the two speeds become equal 
(the intersection points between the broken 
and solid lines), and secondly, this height 
increases with time (= eastward) from about 
2,000 m to more than 5,000 m. This indicates 
that the height of the steering current for this 
typical wave disturbance is dependent on 
time (= eastward longitude). 

The objection might be made that the 
averaging procedure used in determining the 
lower level zonal current values in Fig. 9 
disregards the difference between the air 
motion on the two sides of the frontal zone. 
For this reason the geostrophic wind compo- 
nents on the sea level chart in Fig. 4 b were 
measured. One finds, for the warm sector of 
the Atlantic disturbance 40—so km per hour 
and about the same for the narrow warm sector 
in the disturbance over Germany. Twenty- 
four hours earlier the corresponding warm 
sector values were 40—50 km p.h. over the 
Atlantic and 30—40 km p.h. over Germany. 
Thus the basic result, viz. that the disturbances 
move at a speed considerably in excess of the 


speed of the basic current near the ground, | 


remains valid. 

Near the ground organized vertical motions 
must be of minor importance in as far as the 
horizontal temperature distribution is con- 
cerned, which therefore, in these lowest layers, 
must be controlled primarily by advection. It 
follows that the flat, sinusoidal open warm 
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Fig. 10. Diagram showing the zonal wind speed (com- 

puted geostrophically) over an area surrounding a moving 

wave disturbance. Cy, C2 and C; represent the eastward 

velocities of the disturbance between respectively 03 

GMT and o9 GMT on Nov. 10, 15 and 21 on Nov. 9, 
and 15 and 21 on Nov. 10, 1947. 


sectors of these fast-moving disturbances may 
be considered as a kinematic consequence of 
the fact that the phase velocity of the waves 
exceeds the speed of the basic current near 
the ground (RossBy 1942). In this type of 
‘wave motion the center position at the tip of 
the warm sector will be occupied by ever new 
segments of the front. However, once a small 
occluded frontal portion develops the indi- 
vidual fluid elements which make up the 
occlusion front are unable to move faster than 
the basic current and they must therefore lag 
behind the fast-moving wave center. In the 
formation of this occluded portion of the 
front it may well be that ground friction 
plays a significant rôle, as suggested by 
SCHERHAG (1938). 

Short frontal waves moving faster than the 
zonal wind in the lower troposphere are 
frequently observed. We have included one 
additional example of this, from the second 
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week of November 1947. On the morning of 
November 10 of that year the state of the 
circulation over Europe was rather similar to 
the well-known situation investigated by 
BERGERON and Swosopa (l.c.). A flat, fast- 
moving wave of the “Schnelläufer” type 
passed England during the night between the 
ninth and the tenth with a speed of about 
100 km p.h. For this disturbance a diagram 
(Fig. 10) similar to the one in Fig. 9 was 
constructed, showing the zonal wind speed 
over an area of the same size and shape and 
with the same position relative to the moving 
disturbance as in the case discussed above. In 
Fig. 10 the zonal wind speed is plotted against 
height; the vertical lines marked c,, c, and cg 
represent the speed of the disturbance during 
three six-hour intervals centered around 
consecutive maps. Also in this case the flat, 
open wave disturbance travelled at a far 
greater speed than the zonal wind in the lower 
troposphere. In fact, one must go to the 500 
mb level to find the critical level where the 
zonal wind equals the maximum speed of the 
wave. The wind observations from Liverpool 
for November 9 and 10 (Fig. 11) show an 
interesting distribution of the wind with time 
and height. The wave passed Liverpool about 
midnight and at that time a marked increase 
in the wind at the higher levels was observed. 
The significant changes occurred above sixteen 
thousand feet! Thus the fast-moving open 
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Fig. 11. Diagram showing the distribution with height 
of the observed zonal wind speed over Liverpool for 
November 9 and 10, 1947. 
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Fig. 12. Diagram showing the zonal geostrophic wind 
speed in soo mb for different longitudes during the period 
Feb. 6—12, 1948. The speed is computed with the aid 
of the mean chart in Fig. 3 from the height difference 
between two points, situated 7.5 degrees of latitude on 
each side of the central line of the trajectory belt on Fig. 1, 
one computation being made for every ten degrees of 
longitude. 


frontal waves which frequently have been 
considered as relatively shallow formations, 
are connected with great changes of the wind 
in the upper troposphere, while the corres- 
ponding changes in the lower troposphere 
(apart from the surface layers) are somewhat 
smaller and sometimes hard to detect with 
| the aid of the present wide-meshed aerological 
station net. 

We shall return next to the height varia- 
tions of different isobaric surfaces over Alder- 

rove as recorded in Fig. 5. In view of the 
relatively long time interval between conse- 
cutive soundings it is impossible to fix exactly 
the phase lag between the upper and the lower 
level waves but it may be stated with certainty 
that this lag does not exceed 6 hours, or about 
60°, i.e. the disturbances are almost vertical. 
The temperature variations at the various 
pressure levels have not been entered in the 
diagram as here reproduced, but the data 
show, as might be expected, that temperature 
and height (pressure) vary in phase at all levels 
except at the ground where inspection reveals 
that the correlation between the two curves 
is negative. x“ | 

It was pointed out in connection with the 


analysis of Fig. 9 that the critical level at which 
the zonal speed of the wave disturbance and 
that of the basic current itself become equal 
shifts from about 2,000 m to above 5,000 m. 
In Fig. 12 we have entered the value of the 
zonal geostrophic wind at 500 mb in different 
longitudes as computed from the mean chart 
in Fig. 3. The values calculated refer to a strip 
with a width of fifteen latitude degrees and 
centered around the trajectory belt in Fig.#7: 
This computation indicates quite clearly that 
the basic zonal current itself, at least at the 
soo mb level, decreases markedly eastward. 
The displacement of the critical level upward 
must therefore be interpreted as a tendency on 
the part of individual wave disturbances to 
continue their movement eastward at a some- 
what less decelerated rate than that of the 
basic current in the lower and middle tropo- 
sphere. This result lends some support to the 
impression that the factors controlling the 
behaviour of the wave disturbances are found 
in the upper portions of the troposphere. 

The discussion presented above has brought 
out that the period of the waves is fairly 
independent of longitude and that their phase 
velocity decreases eastward less rapidly than 
the basic current in the lower and middle 
troposphere. With these experimental facts 
as a starting point it is of some interest to 
enquire into the stream line amplitude changes 
that would be observed in a simple transversal 
horizontal wave, superimposed upon a basic 
zonal current the speed of which (U) decreases 
slowly with time. If the transversal wind 
component is given by 


PI y NCOS 


the stream lines may be computed from 


VEIT. ei 
ie nl sin 7 (x — ct) 


where v, is the transversal wind velocity 
amplitude, L the wave length, c the phase 
velocity and y the coordinate northward. 
Since 


24 R. BERGGREN, B. BOLIN AND C.-G. ROSSBY 


where T is the period of the waves, it follows 
that 


The experimental data suggest that cand T 


are relatively constant. If one makes the not 
unreasonable assumption that the kinetic 
energy of the wave motion remains fairly 
constant, the variations in v, must be negligible 
and it follows that the horizontal amplitude 
of the stream lines increases inversely as U 
decreases. 

The rough calculation presented above 
should at any rate serve the purpose of sug- 
gesting that the amplitude increase of occluding 
waves might depend upon a down stream 
decrease in the strength of the basic current 
and that this factor should be explored in the 
final formulation of a theory for frontal waves. 

We shall finally compare the results obtained 
from this study of European frontal waves 
with a somewhat similar study by Hess and 
WAGNER (1948) of certain atmospheric waves 
in the northwestern part of the United States. 
The waves studied by these authors occurred 
with a strong, anticyclonically curved westerly 
to northwesterly current in the middle tropo- 
sphere over the states of Washington, Idaho and 
Montana, changing into a flat and cyclonically 
curved current over the Great Lakes region. 
Station level barograms for a series of points 
from Seattle eastward to Bismarck (North 
Dakota) were plotted for nine consecutive 
days from hourly altimeter settings. The 
resulting diagrams revealed the existence of a 
series of pressure waves which moved with 
great regularity at a speed of about 1 000 miles 
per day (= 67 km per hour) from west to 
east; this speed agrees well with the speed of 
70 km per hour obtained for the particular 
case of wave propagation over western Europe 
studied by us. Aerological data revealed that 
the Pacific disturbances moved. faster than the 
wind in the lower troposphere, the critical 
height sometimes reaching values of up to 
6 km. Radio sonde data analyzed by Hess 
and WAGNER indicated that the disturbances 
west of the Rockies were very nearly in 
phase at all levels, with high positive correla- 
tions between pressure and temperature at all 
heights except near the ground where the tem- 
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Fig. 13. Monthly average speeds of Alberta (SA) and 
North Pacific (Sp) Lows as related to the 10,000 foot 
zonal index (I). (Reproduced from Namias 1947 a.) The 
zonal index is computed for a strip extending over 
North America and limited by 40° and 50°. 


perature variations were insignificant. East of the 
mountains the disturbances moved in over low 
level air masses with a marked north—south 
temperature gradient; the waves then rapidly 
acquired the characteristic low level tempera- 
ture asymmetry of waves on the front 
between polar air masses of Pacific and Con- 
tinental origin. 

In view of the long (12 hours) time interval 
between consecutive soundings no attempt 
was made by Hess and WAGNER to determine 
the period of these waves but an inspection of 
the published barograms for the 6 km level 
suggest a period somewhat in excess of 36 
hours but definitely much less that 48 hours. 

The underlying flow pattern (a strong, 
anticyclonically curved basic current, increasing 
rapidly upward) is very similar in the situations 
investigated by Hess and WAGNER and by us. 
The speed of propagation and the period of 
the disturbances appear to be nearly the same.! 
In both cases a critical level well removed from 
the ground exists and in both cases pressure 
and temperature variations aloft are positively 
correlated. It is reasonable to conclude that 
the perturbations fundamentally are of the 
same character and origin. With these upper 
level similarities are associated marked dissi- 
milarities in the frontal histories of the two 


1 The wave periods obtained from the two series 
studied by Hess and WAGNER and by us should not be 
considered as fundamental parameters. The phase velocity 
of this type of wave motion is intimately connected 
with the speed of the basic current, as may be seen from 
the extremely illuminating curves in Fig. 13. The advec- 
tive motion evidently produces a “Doppler” effect 
which might be eliminated through the introduction of 
a properly defined orbital frequency. 
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Fig. 14. Contours and isotherms on the soo-mb surface, 0300 GMT, February 12, 1948. 


systems. The American disturbances of which 
the Hess— WAGNER series is a particularly 
striking example are often traceable back into 
the Pacific as rapid-moving ‘‘occlusions”. 
After crossing the mountains they acquire 
new fronts and open wave characteristics, 
while the European disturbances start as open 
waves but gradually occlude. It is difficult to 
escape the conclusion that from the dynamic 
point of view the frontal patterns at sea level 
must be regarded as induced phenomena, 
consequences of the fact pd that high level jet 
stream in which the upper waves move, nor- 
mally, but not always, is associated with a 
well-developed front at the ground. This con- 
clusion does not exclude the possibility that 
the upper waves may be intimately connected 


with the powerful, but somewhat more 
diffuse solenoidal field associated with the jet. 
These upper waves would then seem to take 
over röle of the primary pressure waves 
first referred to by von FICKER (1920). A 
somewhat similar point of view has been 
expressed recently also by H. RIEHL (1948). 


II. Blocking and Breakdown of the 
Zonal Motion 


By B. BOLIN 


It was shown in the last section that the 


strong zonal motion which characterized 


the middle and upper parts of the troposphere 
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Fig. 15. Contours and isotherms on the soo-mb surface, 0300 GMT, February 14, 1948. 


over Northwestern Europe and the Atlantic 
during the second week of February 1948 
remained essentially unchanged in spite of the 
occlusion of a series of frontal waves. It is 
the purpose of the present section to describe 
how this apparently stable flow during the 
following weck, between the twelfth and 
twentieth of February, broke down into 
quasi-stationary cyclonic and anticyclonic vor- 
tices. The broad features of this development 
are readily seen from a comparison between the 
synoptic 500 mb charts for February 9, 12, 14, 
16, 18 and 20 (Fig. 4a and 14—18 inclusive). 

On February 9, prior to the breakdown 
(see Fig. 4 a) a well-defined belt of strong zonal 
motion was located in about 55° N and 
extended into Central Europe. Over the 


extreme northern part of Europe the westerly 
winds at 500 mb were weak. South of Europe, 
over the western part of the Mediterranean, 
there was no air motion from the West. Over 
the Baltic this central belt of westerlies split in 
two branches, one extending northward 
towards the White Sea while the other, 
stronger branch flowed southeastward towards 
the Balkans and into Asia Minor. The ultimate 
destruction of the zonal current appears to be 
associated with the gradual intensification and 
slow westward displacement of this split. 
On February 12 (Fig. 14) it occurred further 
to the West, in the region of the North Sea, 
and the process was by then far more striking 
and abrupt than three days earlier. During the 
following eight days the region of branching 
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pushed steadily westward, dividing and de- 
stroying the zonal motion behind it until on 
February 20 (Fig. 18) the process had reached 
Labrador and Newfoundland. By that time 
the jet stream in lat. 55° N over the Atlantic 
had disappeared. In its place one observed a 
series of closed cyclonic and anticyclonic 
vortices, capped to the North by a belt of 
westerlies which extended eastward from 
Greenland at about 70° N, while another, 
apparently continuous zonal belt could be fol- 
lowed from Newfoundland across the Atlan- 
tic through the Mediterranean. It should be 
mentioned that such a zonal wind belt south 
of Europe is observed frequently. There are 
good reasons to believe that this southern jet 
may be identified with the well-defined zonal 


Rh, 


wind maximum at the tropopause level over 
northern India indicated in the normal merid- 
ional cross section for the winter scason 
published in Physikalische Hydrodynamik 
(1933). 

A comparison between the synoptic surface 
charts for February 9 (Fig. 4 b) and February 
18 (Fig. 19) shows that this breakdown of 
the zonal motion aloft was accompanied by a 
corresponding breakdown of the zonal motion 
at sca level. On the latter day the prevailing 
zonal motion over the eastern part of the 
Atlantic and over Europe between about 
35° N and 55° N was in fact from the East. 
This easterly current was surrounded by two 
belts of westerly winds, one to the North of 
Scandinavia and the other over North Africa. 
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The development outlined above exerted 
great influence on weather conditions at the 
ground. The synoptic chart for February 
18 (Fig. 19) provides an excellent illustration 
of the typical weather conditions over Europe 
and the eastern Atlantic following the break- 
down. The most striking feature is the intense 
anticyclone built up over northern Europe. 
To the north of this anticyclone disturbances 
were travelling in the strong westerly current 
and caused rather mild and moist weather in 
Iceland and northern Scandinavia. The tempera- 
ture was between 3 and 6 degrees above nor- 
mal. Most of the precipitation in Scandinavia 
fell on the Western side of the mountain 
range. To the south of the anticyclone cold 
air was brought into Europe from the Arctic 
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Contours and isotherms on the soo-mb surface, 0300 GMT, February 18, 1948. 


Sea and Russia. Thus the temperature decreased 
rapidly between February 14 and 18, and'there 
was a negative temperature anomaly of 5 to 
10 degrees during the rest of the month in 
southern Scandinavia, central Europe, France 
and southern England. Because of the very 
low temperature at the upper levels the air 
was fairly unstable and snow showers were 
frequently observed in this easterly current. 
However, in general the amounts of precipita- 
tion were not large since the moisture content 
of the continental air was low. 

It is evident that the development occuring 
during this period provides a striking illustra- 
tion of blocking as described in recent years 
by Namias (1947 a) and by ELLIOT and SmiTH 
(1948). Namras has presented a particularly 
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Fig. 18. Contours and isotherms on the soo-mb surface, 0300 GMT, February 20, 1948. 


g case of a persistent retrograde 
blocking “wave” in high latitudes (see Fig. 
52 in,Namias 1947). His example, as well as 
most of the cases of blocking studied in the 
last few years have been analyzed with the 
aid of mean charts representing the pressure 
distribution at sea level or at 10000’ for over- 
lapping or consecutive five or seven day 
periods. Even though this method is of 
obvious advantage in bringing out slow 
changes in the large-scale flow pattern it 
also tends to obscure details which may be of 
basic importance. For this reason we shall 
describe below the strongly developed blocking 
wave through its manifestations on the daily 
synoptic charts. It should be remembered, 
however, that individual blocking waves may 
vary considerably in their characteristics. 


illuminating 


Before proceeding to a more detailed 
discussion of the synoptic developments during 
this period we shall attempt to make a nu- 
merical estimate of the westward speed of the 
blocking action. Fig. 14 indicates quite clearly 
that the strong zonal current on February 12 
was limited eastward by a trough line which 
extended from Iceland towards Scotland and 
the North Sea. Similarly, on February 20 
(Fig. 18) the zonal current was limited by 
another trough extending NW--SE and 
located just off the Labrador Coast. From the 
intersections between these two trough lines 
and the latitude circle for 60° N one finds a 
total displacement westward of about 55° of 
longitude in eight days, or an average speed of 
seven longitude degrees per day, in 60° N. 
This rate was not uniform. The maps for the 


m N 


30 R. BERGGREN, B. BOLIN AND C.-G. ROSSBY. 


2 FEB. 18, 1948 
» 0300 GMT 


a 


Fig. 19. Surface map for 0300 GMT, February 18, 1948. 


intermediate days show that the blocking in 
its early stages developed quite rapidly west- 
ward but later on slowed down considerably. 

An accurate estimate from the daily maps 
of the westward speed of the blocking process 
is rendered extremely difficult by the marked 
distortions in the basic soo mb flow pattern 
brought about by a series of rapidly deepen- 
ing waves which during this period moved 
out over the Atlantic from North America. 
The signifidnce of these disturbances will be 
discussed later. To eliminate their effect and 
make possible a determination of the rate of 
retrogression of the blocking process a series 
of overlapping three-day mean contour charts 
for the soo mb level were constructed. Fig. 
20 gives the consecutive positions of a particu- 
lar contour line taken from these charts. One 
finds from this diagram a westward speed of 
four longitude degrees per day in 60° N, or 
about 10 km. p.h. Since the anticyclone which 
developed in the rear of the blocking wave 
changed shape and steadily grew in extent, the 
speed obtained in this manner cannot be 
expected to agree with the first estimate more 
than to the order of magnitude. For com- 
parison it should be mentioned that Namras 
(1947) gives a displacement of about sixty 
longitude degrees per week as a characteristic 
value for high latitude blocking waves. 

It was pointed out above that the advance 
westward of the blocking process did not 
take place at a steady rate but that the entire 
evolution was associated with the rapid 
deepening (and gradual stagnation) over the 


Atlantic of a series of waves which appear to 
have originated in the western part of North 
America. Figure 21 shows the vertical structure 
and approximate period of these waves. The 
data for this diagram were taken from the 
aerological soundings at Caribou, Me (47 N°, 
68° W). The full lines represent, from top to 
bottom, the contour values at 500, 700 and 
850 mb as well as the pressure at sea level, all 
as functions of time. The broken lines give 
the temperature data for the same levels. 
The figure proves that three well-defined 
waves passed the station between February 12 
and February 20 with an average period of 
slightly less that three days. These rapidly 
deepening waves were thus characterized by a 
period almost twice as long as that of the re- 
latively stable waves described in section I of 
this paper. Temperature and pressure varia- 
tions were in phase at the upper levels (i.e. 
the troughs were cold, the ridges warm) but 
at the 850 mb level there was a phase lag be- 
tween the two curves, particularly well-marked 
in the second and deepest wave. In this trough 
the axis of the temperature wave was very 
nearly vertical while the pressure wave at the 
ground was ahead of the upper pressure wave 
by about 24 hours. 

The horizontal amplitude increase of these 
deepening waves is of particular interest in 
connection with the breakdown of the zonal 
circulation. An approximate idea of the 
movements of individual air particles in the 


Fig. 20. Successive positions of the contour line for 5300 

geodynamic meters as taken from overlapping three-day 

mean charts for the soo mb level in February 1948. The 

numbers indicate the periods for which these mean charts 
were constructed. 
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Fig. 21. Contour and temperature (broken lines) varia- 

tions with time at different isobaric levels above Caribou, 

Me. The full lines represent, from top to bottom, the 

contour variations at 500 mb, 700 mb and 8so mb, and 

finally the pressure at sea level. The broken lines represent 
the corresponding temperature variations. 


north—south direction may be obtained by 
studying the changes of the temperature field 
in the soo-mb level. In doing so we assume 
that the particles move along an isobaric 
surface and that the processes are adiabatic. In 
reality these conditions are not fulfilled. How- 
ever, the non-adiabatic processes normally 
produce a cooling of the warm air masses as 
they are carried northward and a warming of 
cold air which is transported southward. 
Furthermore, by and large experience indicates 
that the cold air masses in the lower troposphere 
are subsiding while warm air is rising. Hence 
the effect of vertical movements on the tem- 
perature field at $00 mb is in the same direction 
as that of the non-adiabatic processes and 
both counteract the effect of horizontal (or 
isobaric) advection. It follows that the merid- 
ional displacements must be at least as 
large as those indicated by the isotherm pat- 
terns in the soo-mb level. 

The soo-mb chart for February 12 (Fig. 14) 
shows a narrow, nearly zonal belt around 
so’ N, in which the isotherms are very closely 
packed. Superimposed on this pattern one 
observes a series of fairly evenly spaced iso- 
therm troughs, located (at 50° N) in about 
30° E, 5° W, so’ W and 85° W. The eastern- 
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most of these became important in the final 
stage of the blocking process. For the time 
being we shall consider the development of 
the two western systems. Between Februar 
12 and February 14 (Fig. 15) the isotherm 
amplitudes in these troughs increased mark- 
edly, particularly in the trough located in 
longitude 85° W on February ı2, which 
during these two days moved eaéward about 
forty degrees into the Atlantic just off the 
east coast of Newfoundland. As a result of 
differential motion in the northern and south- 
ern parts of the trough, the horizontal axis 
of the trough took on a marked slant NW— 
SE. During these two days a marked deep- 
ening (pressure fall) occurred in it; this 
deepening in turn gave rise to intense 
southerly winds and warm air advection on 
the eastern side of the trough. Through this 
warm current the stagnant southern portion 
of the next cold trough downstream was 
cut off, to form a cold air pool off the Bay 
of Biscay (Fig. 15). The same intense warm 
current produced pressure rises along its eastern 
boundary which resulted in a broadening to- 
wards the West of the relatively narrow ridge 
of high pressure which on February 12 had 
reached from Central Europe and Scandinavia 
towards the northeastern coast of Greenland. 
On February 14 another deep trough, with 
cold air advection, was located over the Mis- 
sissippi Valley. Two days later (Fig. 16) it had 
reached Newfoundland and Labrador where 
it set off a new and intense air current north- 
ward. The warm air advection and pressure 
rises associated with this current led to a 
complete cutting off of the southern portion 
of the preceeding trough which on Febru- 
ary 16 had been reduced to an isolated cold 
cyclonic vortex off Spain. 

The processes described above were repeated 
a total of four times during the eight day 
eriod under investigation. Each trough, as 
it deepened and slowed down over the central 
part of the Atlantic, gave birth to a strong 
southerly current on its eastern side through 
which the preceeding trough was cut off and 
partially filled. These intermittent warm cur- 
rents gradually broadened the initially quite 
narrow anticyclonic ridge until finally, on 
February 20, a large warm anticyclone 
extended along 60° N from Scandinavia to 
Greenland while the cut-off troughs remained 
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as a series of cold cyclonic vortices on the 
south side of this anticyclone. 

PALMÉN (1949) has stressed the fundamental 
importance of subsidence in the development 
of cut-off cyclonic vortices in the upper 
troposphere. According to him this sub- 
sidence occurs primarily in the northern part 
of the cold air tongues, though there is no 
particular reason why such a distribution of 
the sinking motion should be expected a 
priori. In the series of cut-offs described 
above the required subsidence would be a 
dynamic consequence of the strong, anti- 
cyclonically curved warm current arriving 
from the next trough upstream and would 
occur mainly along the right edge of this 
current. In a large measure, however, the 
initial cutting off scems to result simply from 
horizontal advection. [Compare Staff Mem- 
bers, University of Chicago (1947)]. 

The development discussed above is very 
similar to the seclusion process first described 
by BJERKNES and SOLBERG (1922), though the 
changes in the upper flow pattern studied by 
us are of greater dimensions than the seclusion 
pattern of ordinary short frontal waves. X 

During this period every disturbance trävel- 
ling across the American continent developed 
in conjunction with a wave-like deformation 
of the frontal pattern at the ground. These 
waves were unstable and their amplitudes 
grew rapidly throughout the entire tro- 
posphere. Rapidly occluding frontal distur- 
bances of short wavelength are often observed 
in which no such large-scale changes in the 
upper flow pattern or temperature field are 
found. Section I of this paper contains a 
description of a series of such short waves 
from the preceeding week, just before the 
blocking of the westerlies set in. These short 
waves generated closed cyclonic vortices in 
the lower troposphere below about 600 mb. 
Above this level, however, only minor waves 
were observed and the individual air particles 
performed an oscillating motion with an 
almost constant amplitude while travelling 
eastward. No closed vortices were formed by 
the cutting-off of cold air in the manner 
described above. Consequently these shorter 
waves were not unstable in the sense of producing 
a complete breakdown of the wavelike flow pattern 
at the upper levels. It is then necessary to observe 
the difference between these two kinds of 


waves in the atmosphere, both of which appear 
as occluding frontal waves at the ground but 
represent quite different processes aloft: 1) short 
frontal waves (wave length 2 500 km or 
less) which may occlude but nevertheless do 
not change the upper flow pattern and 2) some- 
what longer, apparently unstable waves which 
develop in conjunction with a complete 
breakdown of the zonal flow. 

The reservation implied in the term “appar- 
ently unstable” is made advisedly. A generally 
accepted theory for the instability of these 
somewhat longer waves has not been developed 
yet. Assuming however, that such a theory 
could be formulated it would presumably 
lead to the conclusion that a certain range of 
wavelengths is unstable; it might even be 
possible to show that there exists a particular 
wave length of maximum growth rate, as was 
done by Eapy (1948) for the particular atmo- 
spheric model studied by him. By dropping 
the restriction to small amplitudes it might 
be possible to show that the speed of propaga- 
tion of the individual troughs would slow 
down as the amplitudes grow, and that adja- 
cent waves would influence each other. How- 
ever, such a theory could not by itself account 
for the characteristically uneven spacing of 
the troughs that developed during the period 
studied by us and that expressed itself as a 
tendency of the troughs to pile up off the 
coasts of Europe. Any theoretical treatment 
of the growth and breakdown of a series ‘Of 
small, sinusoidal perturbations would inévi- 
tably lead to identical behaviour of the indi- 
vidual troughs and would not by itself be 
able to account for the marked shortening in 
wave length observed over the eastern part of 
the Atlantic. One is forced to conclude that 
this observed decrease in the spacing between 
consecutive troughs as they come in from the 
Atlantic, inexplicable on the basis of instability 
in the commonly accepted sense, must be 
looked upon as a manifestation of the blocking 
process, whatever its cause, and must be 
attributed to the successive (rather than si- 
multaneous!) retardation of each consecutive 
trough in a given region. Under those circum- 
stances it 1s impossible to exclude completely 
the possibility that the observed amplitude 
increase of the horizontal isotherm pattern in 
each trough (ie. the instability) somehow 
may be the consequence of the retardation 
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of the troughs, rather than the cause. One is 
once more led to the same conclusion as in 
the previous section, viz. that the theory for 
these different types of perturbations should 
take into consideration the variation with 
longitude of the underlying basic current 
pattern itself. 

In spite of the complications pointed out 
above it may be of some interest to explore 
the applicability of existing theories for un- 
stable waves to the deepening disturbances 
analyzed in this section. The first studies of the 
dynamics of long waves were carried through 
for an incompressible and barotropic atmo- 
sphere to clarify the röle played by the varia- 
tion of the Coriolis’ parameter with latitude 
in determining the scale of atmospheric 
disturbances. In view of the restrictive assump- 
tions on which these studies were based they 
could not provide any information concerning 
the stability of atmospheric waves. Subsequent 
investigators attempted to solve the stability 
problem by taking the barocline state of the 
atmosphere into account. JAW (1946) studied 
the problem in two dimensions by introducing 
the effect of the horizontal solenoidal field 
and found a critical wave length beyond which 
the waves become unstable. A more general 
attack was made by CHARNEY (1947) who 
treated the three-dimensional problem and 
took the vertical shear into consideration. He 
derived an instability criterion which indicated 
that the instability increases with vertical 
shear, vertical temperature lapse rate and 
latitude but decreases with wavelength. 

When applying these theories to the upper 
air charts it is necessary to keep in mind that 
they .were derived for a rather simplified 
atmosphere and by assuming small pertuba- 
tions, in order to permit linearization of the 
hydrodynamic differential equations. Thus, 
while the theories may give the conditions 
necessary for the initiation of instability they 
cannot describe in what manner the subsequent 
development takes place. In spite of this, it is 
of interest to compare the initial conditions 
in February 1948 with those prescribed by 
the theory. The mean value of the vertical 
shear has been computed for each day in a 
1 soo km broad strip along the jet between 
longitude 130° W and 60° W, the region 
where the waves grow most rapidly. The 
figures thus obtained (table II) show an in- 


crease of the shear as the development of the 
blocking starts, indicating a qualitative con- 
firmation that the waves have become un- 
stable. 


Table I. Mean vertical shear in a 1500 km 

broad band along the jet between longitude 

130° W and 60° W, computed from the isotherms 
in the daily 500-mb charts for 0300 GMT 


Mean vertical shear 


pate m/sec/km 
ÉCDE AT ne 1.7 
12 Faro 00 Oto rar I.9 
ee ER 2.2 
LAS Rene chlore 2.3 
Re EE Tee 24 


The magnitude of the shear, however, is 
considerably greater than is necessary in order 
to get instability according to CHARNEY. The 
computations represent a mean value for a 
rather broad band. Nevertheless, the actual 
shear so determined is almost twice as great as 
the critical value demanded by Charney’s 
theory: 1.0 m/sec/km. In the frontal zone, 
where the strongest concentration of iso- 
therms exists, the shear is still greater: 4—6 
m/sec/km. Charney himself reached the con- 
clusion that the westerlies are the seat of a 
nearly permanent instability. Since a complete 
breakdown of the waves in the manner 
described by us occurs fairly infrequently it 
appears impossible to explain this development 
as the consequence of a nearly permanent 
instability of the zonal flow. Once more the 
impression is gained that the major change in 
the flow pattern associated with the blocking 
process may be the primary factor and that 
the amplitude of the individual waves grow 
as a result of the blocking and of the retarda- 
tion in phase velocity produced by it. In this 
connection it is instructive to consider the 
changes in the longitudinal distribution in 
zonal motion at soo mb which are observed 
during the period (Fig. 22). In this figure we 
have plotted the zonal motion at soo mb be- 
tween 45° N and 65° N as a function of longi- 
tude for three slightly overlapping three-day 
periods. The westward progress of the blocking 
is brought out by the diagram as is the fact 
that the zonal motion throughout the entire 
period decreases in strength towards the east. 
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Fig. 22. Distribution at 500 mb of the zonal geostrophic 

motion between 45 and 65° N as a function of longitude, 

at different stages of the blocking process. The numbers 

indicate the partly overlapping three-day periods for 
which the computations were made. 


At least superficially the sharp break in the 
zonal current at the western boundary of the 
particular blocking wave described in this 
paper resembles the standing wave or hydraulic 
jump which develops in open channel flow 
when the depth of the water current drops 
below a certain critical depth. The existence 
of a critical relationship between speed and 
depth may be derived from the laws of con- 
servation of mass and momentum but can not 
be established by means of linearized equa- 
tions. Fig. 23, reproduced from ADDISON 
(1944) gives an idealized representation of 
the flow through such a standing wave. One 
finds that the total energy is less at section 2 
than at section 1; the surplus is transformed into 
violent turbulent motion in the standing wave. 
If this analogy is accepted the large-scale 
turbulence which appears in the present case 
as the waves break down would be the con- 
sequence of the formation of a blocking wave. 
The principal problem is then to determine 
the conditions under which such a change in 
the basic current takes place. In an investiga- 
tion of this process it may well be necessary 
to take into consideration the three-dimen- 
siomal structure of the jet stream even though 
it has been shown (YEH 1949) that a consider- 
able amount of information concerning the 
blocking process may be obtained from the 
analysis of barotropic models. 

We shall now consider the final stage in the 
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development of the high latitude blocking 
anticyclone. It has been shown above that the 
gradual westward growth of this system was 
associated with a series of intermittent warm 
southerly currents emanating from successively 
deepening Atlantic troughs. Towards the end 
of the period (February 17—20) this anti- 
cyclone was .cut off completely from the 
warmer air in southerly latitudes by a major 
change in the flow pattern over Finland and 
Russia. As a result of the strong pressure rises 
over the northernmost part of the Atlantic 
the contours to the east of the ridge took on 
an almost straight meridional orientation. The 
cold air over the Arctic Ocean then had the 
opportunity to move southward. In the eastern 
part of this current the cold air advection 
gave rise to rapid cyclogenesis, but a portion 
of the cold air was deflected westward and 
appeared over northern Germany and southern 
England as a very intense easterly jet. With 
this development the breakdown of the 
westerly flow had been completed. On Feb- 
ruary 20 (Fig. 18) the zonal index at the 500 
mb level, computed for the region between 
45° and 65° N and between 60° W and 40° E, 
had dropped to — 3.5 mps. To get an idea 
of the structure of the warm anticyclone and 
of the cold vortex on its south side a cross 
section has been constructed extending from 
Jan Mayen in 70°N to Maison Blanche 
(Algiers) in North Africa (Fig. 24). From this 
section, as well as from the individual $oo-mb 
charts one gets a strong impression of the 
large-scale horizontal mixing associated with 
the breakdown process, warm tropical air 
having been brought far to the north while 
cold arctic air has been drawn into the temper- 
ate and subtropical zones. 
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Fig. 23. Section through a standing wave (hydraulic 
jump) in a water current. (Compare e.g. ADDISON, 1944.) 
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Fig. 24. Vertical cross section from Jan Mayen to Algiers (Maison Blanche) for February 20, 1948, 0300 
GMT. Frontal boundaries, inversions and tropopause are indicated by thick solid lines when distinct, and by 
thick broken lines when indistinct. Thin solid lines represent isentropes and thin broken lines isotherms. 


Some interesting features of the easterly jet 
may be pointed out. The relative vorticity of 
an individual air particle in the current (at 
500 mb) decreases in spite of the southward 
movement. The complete vorticity equation 
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Fig. 25. Distribution of wind speed in the easterly jet 

over southern England on February 20, 1948, 0300 GMT. 

Thin solid lines indicate the wind speed normal to the 

section and thin broken lines the isentropes for dry air. 

Frontal boundaries, inversions and tropopauses are 
given by the thick lines. 


shows that this decrease may be explained as 
an effect of horizontal divergence within the 
current, of a horizontal solenoidal field or of 
certain vertical movements. In the present 
case the solenoids act in the opposite direction 
to the one required, i.e. they tend to produce 
cyclonic vorticity, as may be seen from the 
indicated cold air advection. There are two 
observations which suggest that the decrease 
in vorticity was brought about by divergence. 
In the first place, the height of the tropopause 
over the cold current dropped to lower and 
lower levels, reaching a minimum of 4 500 m 
over Paris on February 20 (see Fig. 24). In 
the second place, from an approximate analysis 
of the vertical motions it can be shown that 
there was no general descending motion 
in the lower troposphere, below about 2 000 
m. The last statement is borne out also by 
low-level convection in the cold air referred 
to earlier. One is thus led to the conclusion 
that the upper portions of this cold air current 
must have spread out laterally during the 
movement southward. 

The wind distribution in the jet is shown 
in Fig. 25 which is based on observed winds 
extrapolated by means of the temperature 
field and the geostrophic assumption. On 
February 20 the jet was located over southern 
England with a pronounced maximum of 
about 70 mps in the 400 mb level. The maxi- 
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Fig. 26. Idealized sketches 


of the development of unstable waves at the soo mb level, in association 


with the establishment of a blocking anticyclone in high latitudes. Cold air in blue colour, warm air in red. 
Solid lines are stream lines and broken lines the frontal boundaries. 


mum anticyclonic shear which is likely to 
be observed in a zonal current, except as a 
transient phenomenon, should correspond to 
constant absolute angular momentum. In this 
case both observed and computed winds seem 
„to indicate that this critical value was reached 
on the northern side of the easterly jet. In this 
connection it should be mentioned that pro- 
fessor J. BJERKNES in an as yet unpublished 
investigation has attributed :the anticyclonic 
sweepback of deep upper troughs of the type 
here studied to the presence of inertia in- 
stability on the south side of a westerly jet 
as it crosses the upstream ridge. 

There are two additional factors which 
appear to be of importance in connection 
with the establishment of the easterly jet. 
As a result of the blocking process one branch 


of the zonal westerlies had been shifted to 
high latitudes. The stationary wave length of 
planetary waves, measured in longitude de- 
grees, is given by 
À = 180 Vos , 
a Q cos ¢ 

where U is the velocity of the zonal current, 
a the radius, © the angular velocity of the 
earth (in c.g.s. units) and @ the latitude. A is 
the wave length measured in degrees longi- 
tude. In the derivation of this formula the cur- 
vature of the earth was not taken into account 
and the expression therefore can not give more 
than an approximate value for the wavelength. 
Nevertheless, it is evident that the stationary 
wavelength is greater in high than in low 
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latitudes for a given value of the zonal current. 
This difference in wave length between high 
and low latitudes shows up rather clearly on 
the soo mb charts for the latter part of the 
period. It appears safe to conclude that the 
ultimate establishment of a quasi-stationary 
cold trough over the White Sea region may 
be regarded as the consequence of the long 
quasi-stationary wave length in high latitudes 
and of the development of a strong southerly 
current in the region of Greenland. Through 
this quasi-stationary White Sea trough cold 
arctic air was brought far to the south. It was 
shown earlier that the westward deflection 
of this cold current was associated with marked 
horizontal divergence. However, it is not 
unlikely that the marked anticyclogenesis on 
the south side of the developing northern 
zonal current may have been a significant 
contributing factor, establishing air motion 
from the east over central Europe prior to 
the arrival of the cold arctic air. This appears 
particularly likely from the character of the 
flow pattern at the soo mb level on February 18. 

NAMIAS (1947 b) has described the establish- 


ment of a very intense blocking anticyclone 
in February 1947. This period seems also to 
have been characterized by a development 
similar to the one described in this paper. 

In order to sum up the various essential 
features of the blocking process described 
above we have prepared the idealized sketch 
maps reproduced in Figure 26. In these 
sketches stream lines are solid and the dash- 
ed lines represent the frontal zone. Cold air 
is indicated by blue colour, warm air by 
red. The width of the frontal zone on the 
sketches is inversely proportional to the in- 
tensity of the front. In reality the front is not 
as sharp as shown in these figures but this 
overemphasis was made to give a clear picture 
of the principles involved. For the same 
reason the anticyclones to the north have been 
drawn as separate vortices though they actually 
appear to merge into a single system. The 
retrogression, the cutting off of cold cyclonic 
vortices and of warm anticyclones as well as 
the final injection of a cold current from the 
northeast have been incorporated in these 
figures. 
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Abstract 


It is shown that the large-scale quasi-stationary disturbances of the middle-latitude westerlies 
are produced by the forced ascent of the westerly current over the continental land masses. 
Friction is found to have an important modifying effect on the motion. Using these re- 
sults a numerical method is devised for predicting the height profile of the soo mb pressure 
surface at a fixed latitude. The method involves the use of the notion of an “equivalent- 
barotropic atmosphere” and of the geostrophic approximation. Six actual forecasts are made 
for a period of one day and the results compared with observation. The accuracy obtained 
is thought to justify incorporation of the method into day-to-day forecast procedures. 


I. Introduction 


In an article by one of the co-authors,? a 
program for numerical weather prediction was 
outlined in which it was proposed to consider 
a hierarchy of atmospheric models whose 
study would lead to an increasing com- 
prehension of the physical and numerical 
aspects of the forecast problem. The most 
elementary model was a barotropic atmosphere 
in which the motion is regarded as consisting 
of small perturbations on a zonal current. The 
problem of forecasting these perturbations 
constitutes the simplest non-trivial instance 
of a numerical forecast problem. It is the 
purpose of the present article to discuss this 
case as a step towards the realization of the 
general program. It is also hoped that the treat- 
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ment of motions corresponding closer to 
reality than those previously dealt with will 
prove of value to the forecaster in the field. 
In previous investigations of this barotropic 
model, kinematic constraints were always 
imposed in order to permit the use of simple 
mathematical techniques. For example, the 
infinite wave train on a plane earth served as 
a model for the upper air cyclone trough and 
the point source as a model for the localized 
action of solenoidal fields. C.-G. Rossy 
(1945), pointing out that the study of purely 
periodic or isolated systems does not suffice to 
explain phenomena associated with energy 
transfer, introduced the concept of dispersion 
into meteorology. But this concept too is 
applicable only to a restricted class of mo- 
tions, namely motions that can be described 
as wave trains with slowly varying wave- 
lengths and frequencies. In order to deal with 
the most general types of interaction between 
systems, it is clear that all kinematic constraints 
must be dropped. In so doing, numerical 
methods must replace the purely analytic 
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mathematical techniques, which no longer 
suffice to deal with the complexities of the 
motion. Although there is a loss in elegance 
of treatment in using numerical methods, the 
gain in generality more than compensates; it 
becomes possible — for the first time — to 
treat actual initial and boundary condition.s 

A preliminary step was taken in (NWP) 
where a method was given for predicting the 
actual height profile of the soo mb pressure 
surface at a prescribed latitude on the assump- 
tion that the disturbances could be regarded 
as perturbations of infinite lateral extent on 
a uniform zonal current. The agreement of 
an actual 24 hour forecast with observation 
was better than had been expected, but there 
was also a marked discrepancy. A deep and 
extensive trough in the western Pacific and a 
wedge in the eastern Pacific were predicted 
to move rapidly westward but were observed 
not to move at all. 

The persistence of the Pacific trough and 
wedge system suggested the explanation that 
tne motion is composed of a traveling “free” 
perturbation and a permanent “forced” per- 


turbation produced by the action of geo- 
graphically fixed perturbing forces. The 


authors were thus led to a study of the 
quasi-permanent atmospheric perturbations. 
It was found that the observed mean dis- 
turbances at the soo mb level could be ad- 
equately explained as an effect of the forcing 
of a zonal current over the continental eleva- 
tions, providing surface friction and the lateral 
variation of the motion were also taken into 
account. 

Having been convinced that the main dis- 
crepancies in the forecast could be elimin- 
ated by taking into account topographical 
and frictional factors, the authors devised 
a method for incorporating these factors into 
the forecast equations. Applying the method, 
24 hour forecasts were prepared for the period 
Jan. 8—13, 1946, and were compared with 
observation (see Figs. 8—13). In most cases the 
accuracy appears good enough to justify in- 
corporating the method into routine fore- 
casting, and in no case are the errors as great 
as in the original forecast (Fig. r). The method 
is objective and can therefore be used to supple- 
ment forecasts of a more subjective character. 
In certain cases it is believed to yield a better 
forecast than the standard method of extra- 


polation, particularly where deepening or 
filling due to horizontal energy dispersion 
takes place. 


II. The equivalent-barotropic atmosphere 


It is a matter of experience that the large- 
scale perturbations in the atmosphere have the 
character of “external” waves in which the 
shape of the streamlines (or isobars) is approx- 
imately the same at all levels, and where the 
increase of wind with height is similar along 
all verticals. These properties make it possible 
to deal with the large-scale motion approx- 
imately as a two-dimensional problem by 
using a method given in (NWP). The deriva- 
tion of the method will be indicated briefly 
here. 

The vorticity equation can be written with 
good approximation 


0 
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in which the operator 7, denotes the indivi- 


dual rate of change with time. & is the ver- 
tical vorticity component relative to the carth, 
f is the Coriolis parameter, u is the compo- 
nent of velocity in the x-direction (cast), v is 
the component of velocity in the y-direction 
(north), @ is the density, Q is the angular 
velocity of the earth and @ the latitude. Here 
€ is considered as small compared to f, and 
the solenoid term is combined with the diverg- 
ence term to give the mass-divergence. Form- 
ing the average (denoted by a bar) of the 
vorticity equation in the vertical direction 
with respect to pressure and utilizing the ten- 
dency equation, we find 


ae 
A , 


where the subscript ‘‘o” is used to denote 
surface values, and H = RTy/g is the height 
of the homogeneous atmosphere. w is the 
vertical velocity component, p is the pressure, 
T the temperature, R the specific gas constant, 
and g the acceleration of gravity. It should 
be noticed that the effect of vertical solenoids 
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does not occur explicitly in this equation, 
which, so far, involves only a slight approxima- 
cion. It is now assumed that the velocity field, 
with sufficient approximation, can be written 


u=u(x 7, 9 A(p)v =v (yA (p) (3) 


(hence A = 1), and that the vertical motion 
for the large-scale systems is so small that its 
effect can be neglected in the individual de- 
rivative of the vorticity. With these assump- 
tions, (2) can be rewritten as 
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According to the assumption (3), there must 
be a certain level p where Alp) =A =1; at 


this level we have # =u, v =v, and € =€, 
so that equation (4) may be written 
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at the level p = p. In a barotropic fluid, we 
have A(p)= 1, so that equation (s) (with 
A? = 1) holds at all levels. Equation (5) can 
therefore be interpreted as stating that an 
equivalent-barotropic atmosphere can be defined, 
whose motion corresponds to the motion 
of the baroclinic atmosphere at a certain 
level. This level (p = p), where equation (5) 
en will be called the equivalent-barotropic 
evel. 

According to CHARNEY (NWP) a proper 


value of A? in the atmosphere is 5/4 and the 


value of p lies between soo mb and 600 mb; 
the value soo mb will be chosen here as a 
matter of convenience since the data at this 
level is given explicitly in synoptic reports. 


III. The quasi-geostrophic approximation 


The large-scale perturbations in the atmo- 
sphere show a rather slow speed of propaga- 
tion relative to the air. This fact indicates that 
these perburbations are of the type called 


“planetary” by RossBy. Besides such planetary 
motions, however, the equations of motion 
permit solutions which may be characterized 
as combined gravitational, inertial, and sound 
waves; such waves will have a speed of pro- 
pagation comparable with the speed of sound.! 
The fact that the equations permit as solutions 
such fast-moving perturbations will cause 
certain difficulties for numerical integration of 
the equations. CHARNEY (1948) has shown, 
however, that by making use of the geo- 
strophic approximation one can derive a 
system of equations where the fast-moving 
perturbations in which we are not interested 
are eliminated. The geostrophic approxima- 
tion is applied after having eliminated the 
horizontal divergence in the vorticity equa- 
tion by means of the equation of continuity. 
ELIASSEN (1949) derived virtually the same 
equations by utilizing the geostrophic approx- 
imation in a somewhat different manner. Using 
Cuarney’s method, we simply introduce the 
geostrophic approximation into the vorticity 
equation (5) for the equivalent-barotropic 
level. This gives 
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where B = df/dy, ug, vg are the components 
of the geostrophic wind, and ¢, the geostrophic 
vorticity, at the equivalent-barotropic level. 

In the following, we shall deal with the 
motion in an isobaric surface and shall therefore 
replace the pressure as dependent variable by 
the height z of the isobaric surface. We 
then have 
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The geostrophic vorticity is, 
approximation, 
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1 It should be emphasized that we are here concerned 
with very long, external waves only, and that such small- 
er scale phenomena as frontal waves are disregarded. 


with good 
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The pressure tendency is 


ry CET: (9) 


When (7), (8), and (9) are substituted into the 
vorticity equation (6), an equation in z alone 
results, except that it contains the tendency 
and vertical velocity at the ground on the 
right-hand side. In the case of motion over 
level ground, wy =0, but the tendency at 
the ground still remains. When the computa- 
tions were started, this indeterminacy was 
removed by making the assumption = = = 
¢ 
It was later recognized that it is more con- 
sistent with the assumptions made previously 


402 Oz 3 ; 
to write Er DO where x is the ratio of 
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the surface geostrophic wind to the wind at 
the equivalent-barotropic level. However, the 
difference is insignificant, since the whole 
term may be dropped without noticeable 
effect on the results. 


IV. A preliminary forecast attempt 


The motion at the equivalent-barotropic 
level (soo mb) is considered as consisting of 
small perburbations superimposed upon a 
constant zonal current U. For the time being, 
the height perturbation z of the soo mb 
surface will be considered as independent 
of y. Assuming motion over a level ground 
(w, = 0), the quasi-geostrophic vorticity equa- 
tion (6) for the equivalent-barotropic level 
becomes 


B à ge) oe Le pe 


— 2-0: NEO 
ot dx) dx2 | Ox ot ( ) 


where 


Mk 


A? on (11) 


Strictly speaking, the factor U’ = AU =5/4 U 
should be used in place of U in (10). How- 
ever, it was not found to be worth while to 
make this distinction, since U is not uniquely 
determined from the data in any case, but 
depends on the method of averaging. 
Equation (10) is assumed to hold at a re- 


presentative latitude. It is convenient to use I 
day as the time unit and the radius of the 


latitude circle as the unit of length. The 
solution z (x, f) must then be considered as 
periodic in x with the period 27. By means 
of FOURIER analysis we find the solution in 
terms of the initial distribution z (x, 0) in 
the form 

zZ UE st) = 


== 2 (x, 0) + fe (a, 0) I» (x — ax, t) da (12) 


where the influence function (or GREEN’S 
function) I (x, t) is defined by the FOURIER 
series 


I Es = inbt . 
Le Ge t) = Le > ee Ee 1) ee. (13) 


n= — © 
in which 


b=B+22U (14) 


In the units used (x in radians of longitude, t 


in days) we find 
B = 4x cos? y, 


0272 


DE sin? 2p = 2.5 sin? 2, (15) 


where r denotes the earth’s radius. 
Taking 9 = 45° as the representative latitude, 


(16) 


The factor b varies with the current velocity, 
but this variation is very slight since U is 
always of the order of 1/3 radian per day. 
One may therefore consider b as a constant. 
In the computations the value b =7 was 
adopted. 

The series (13) for the influence function 
converges rather slowly, so that the straight- 
forward computation by evaluating a sufficient 
number of terms in the series would be a very 
laborious task. However, the authors’ associate, 
G. Hunt, has given a method by which the 
computation can be done without too much 
work. Equation (13) may be written 


br Mes 
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where 


übt 
LED ES (em —x) em (8) 


n=—o 


is the influence function for non-divergent 
planetary motions. The first sum on the right- 
hand side of (17) converges so rapidly that 
it is sufficiently accurate for our purpose to 
sum from — 20 to + 20, as an estimate of the 
remainder shows. As to the function Jy, (18) 
may be written 


1 + co 4 ( )2 ICO Ds 
ibt eu ibt ee 
LEE 27 > À 27-2! > n? 
Nn=—o N= — © 
+» L + © j 
4 (ibt)? eine (ibt)4 eine ei 
27-3! nr 27+ 4! n* 
a ibt s 
I a ibt (ibt)? 
Gi 27 E : É an? F 
N—= — oo 
(CD ACL ES 
nr 3 4!n* a (19) 


The four first series can be summed exactly; 
they represent polynomials in x of the first, 
second, third, and fourth degree, respectively. 
The number of such terms to be taken is a 
question of economy and depends on the 
accuracy wanted. The convergence of the 
last series is so rapid that summation from 
— 20 to + 20 will give more than sufficient 
accuracy. 

In this way, the functions I,(x, t) and Ie(x ft), 
were computed for £ = 1 day accurately to 
two decimal places.“ The results are given in 
table I, and I,,(x, 1) is graphed in Fig. 5. 
The prognostic formula (12) was tried out 
on an example from the carefully pre- 
pared and well-analysed Historical Weather 
Maps series published by the Air Weather 
Service of the U.S. Air Force. The height 


1 A more accurate computation of a GREEN’s func- 
tion for non-divergent planetary waves, for several diff- 
erent latitudes and time intervals, has been made inde- 
pendently by G. E. FORSYTHE in a paper to be published 
in the Journal of Research of the National Bureau of 
Standards. 
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Fig. 1. The upper curve gives the observed height pro- 

file of the soo mb surface at 45° N for January 12, 

1946, 0400 GMT. The lower continuous curve gives 

the observed change in the following 24 hours, and the 
lower dashed curve the predicted change. 


of the soo mb surface along 45° latitude for 
January 12, 1946, 0400 GMT was taken as 
the initial distribution of z; and the correspond- 
ing profile for January 13, 1946, 0400 GMT, 
was computed by means of equation (12) 
with U = 20° longitude per day. The integral 
in (12) was approximated by summing over 
intervals of 10° longitude. The profile of the 
soo mb surface along 45° latitude on January 
2, and the computed and observed 24 hour 
changes are shown graphically in Fig. 1. 
The zonal profile of the s00 mb surface may 
be described as an almost stationary, very long 
wave pattern, upon which is superimposed a 
system of shorter, migrating waves. The long 
wave, stationary pattern consists of two main 
ridges at the west coasts of Europe and North 
America, and of troughs over the North 
American east coast and the western Pacific. 
The computed forecast is not bad, as far as 
the shorter waves are concerned; but the 
formula gives a rapid propagation of the 
very long wave pattern toward the west, in 
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Table I 
ee ee ee See EEE AE 
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(degrees [Io(x, 1) Io(-x,1)|Iz.5(x,1) I2.s(-x,1)|[Is(x,1) Is(-x,1) To(x,1) Iuo(-%,1)|T1a(&,L) Hal-%,1)|Is(x, 1) Iye(-x,1) 


long.) 


mn nn nn m m Mn en et 2 ea a a FE 


[6) 20. 20,30 | — 6.64 0.36 — 5.78 
Io = a 0.34 | — 3.44 0.44 2.05 
20 — ere O53 ot |e— 106 0.50 LL 
30 O.14 0.43 | — 0.68 Oss —OV4T 
40 0.78 0.47 | — 0.17 0.54 = Où 
50 0.38 0.50 0.06 0.54 — 0.02 
60 0.94 0.52 0.16 0.53 0.02 
70 0.74 0.53 0.18 0.51 0.03 
80 0.480 0353 0.17 0.48 0.03 
90 0,22 20.53 0.16 0.44% 0.02 

100 — Oe | 0.48 0.14 0.41 0.0I 
TIO Oa. 0.43 0.12 (PR 0.00 
120 sO. Om 10.27 O.rt 0.34 0.00 
130 0745) 20.28 0.10 0.30 OO 
140 O6 Ole TO O.TE 0.2 — OL 
150 =— 9.47 0.09 OF 12 0.2 CO 
160 — 0,42 —0.02 0.13 0.21 — 0.01 
170 — 0/34 0:74 O.14 O.18 0.00 
180 — 0.24 — 0.2 0.16 0.16 0.00 


99999999999999909HH 


22 |—5-32 1.68 |—5.06 1.94 |—4.88 2.12 
Io |— 2.17 {35 11:86 1.35 = 71,64 E36 
98 |— 0.86 1.02 |—0.69 0.96 |—0.55 0.88 
86 | — 0.30 0.76 |— 0.22 0.64 |—0.16 0.53 
72 — 0.07 0.56 — 0,04 0.44 = 0.02 0.33 
58 |—o.o1r O.41 — O.01 0.28 0.00 0.20 
47 0.02 0.29 0.01 0.19 0.02 0.12 
38 0.01 O.21 0.00 ORT 0.01 0.07 
30 0.01 0.15 0.00 0.08 0.00 0.04 
2 0.00 0.10 020% 0.05 0.00 0.02 
18 0.00 0.07 — 0.01 0.04 0.00 0.01 
14 |— 0.01 0.05 — 0.01 0.02 0.00 0.00 
ko —K SRG 0.03 D OO 0.07 0.00 0.00 
OFF 220,0 F 0.02 0.00 0.00 0.00 0.00 
05 | — 0.01 0.02 0.00 0.00 0.00 0.00 
03 0.00 0.01 0.00 0.00 0.00 0.00 
OI 0.00 O.01 0.00 0.00 U,00 0.00 
00 0.00 0.00 0.00 0.00 0.00 0.00 
00 0.00 0.00 0.00 0.00 0.00 0.00 


The influence function Ia® (x, 1) for a = 0, 2.5, 6, 10, 14, and 18. 


contradiction to the observed persistence of 
this system. . 

As an explanation for the discrepancy, it is 
assumed that the very long waves are forced 
perturbations, set up by some geographically 
fixed disturbing forces. There are two conceiv- 
able ways in which the motion can be effected 
by the underlying surface: by thermal action, 
or by the forcing action of continental eleva- 
tions. The former has been suggested by 
Rosssy (1939) and Haurwirz (1940 Il), but 
is here tentatively discounted on the groun 
that the reversal in phase of the fixed perturba- 
tions, which should be expected to accompany 
the reversal of the heat and cold sources 
from winter to summer, is not observed to 
occur. The latter is therefore assumed to the 
dominating effect, and we are led to consider 
the nature of topographically forced perturba- 
tions. 


V. The effect of continental elevations 


Instead of putting the vertical velocity at 
the ground in equation (6) equal to zero, we 


oh oh 


now write wy = Un — + Ya, Where h is 
Ox oy 


the height of the ground above mean sea level. 


The perturbation equation for stationary 
motion independent of y becomes 


d3z dz dh 
UT | Lone et dx’ 


(20) 


where h(x) is the mountain profile along the 
latitude circle in question. Assuming Up, the 
geostrophic zonal current at the ground, to 
be a certain fraction x of the velocity U at 


Table II 
An 
N 

@=10 ee eier | a? =18 
3 10-070 OO | 0,087 
2 — 0.100 — 0.080 —0.064 
="; — 0.506 — 0,433 —0.382 
oO — 0.212 — 0.182 —0.161 
I G.305 0.314 0.326 
2 O.119 O.112 0.102 
3 O7 O.149 0.123 
4 0.065 0.051 0.038 
5 0.095 0.065 0.047 
6 0.034 0.022 0.014 
7 0.049 0.030 0.016 

8 0.017 

9 0.023 


Numerical constants A» in the forecast equation (45). 


44 J. G. CHARNEY AND A. ELIASSEN 


the soo mb level, we get after dividing through 
by U and integrating: 


C2 Wie : 
al — x}2h(x) (21) 


From this equation it follows that over the 
oceans, where h = 0, the stationary perturba- 
tion must take the form of a pure sine wave, the 


wave length of which is 2m VU/B, RossBy’s 
stationary wave length. 

This result may be compared with the nor- 
mal pressure profile at 20,000 feet along 45° 
latitude for the month of January, taken from 
the Normal Weather Maps series of the U. S. 
Weather Bureau. This normal profile, which 
is shown converted to height units at soo mb 
in Fig. 4 (solid line), may be considered as 
approximating the stationary profile. As 
expected, the profile is roughly sinusoidal 
over the oceans; but the wave length is 140°— 
150° of longitude, whereas an estimate of the 
zonal wind speed shows that Rosspy’s sta- 
tionary wave length is about 80°. We are 
therefore forced to the conclusion that (21) 
is not the correct mathematical expression of 
the stationary motion. 

A more correct explanation of the stationary 
pattern follows from the consideration that 
the lateral variation of a topographically 
forced perturbation must be determined by 
that of the continental elevations, and these 
are more closely approximated in meridional 
section by the crest of a sine wave than by a 
horizontal line. Accordingly we allow for 
the meridional variation by assuming that the 
perturbation has a sine variation in the north— 
south direction with a half-wave length cor- 
responding roughly to the effective lateral 
extent of the land profile across the zonal 
current, about 25° of latitude. The motion is 
therefore of the type studied by Haurwirz 
(1940, I). Putting 


in the expression for the relative vorticity (8), 
we obtain, instead of equation (21) 


’ 


02z 


Ox? 


sz = — xA*h(x) 


where 


st = m? (24) 


In the length unit used (degrees of longitude 
at 45°) s represents the number of stationary 
waves encircling the 45° latitude circle and 
m = n/D, where D is the north—south half 
wave-length of the perturbation. The value 
of s must be approximately 2.5 to correspond 
to the stationary wave length of 140° longitude 
observed over the oceans (see Fig. 4). The 
value of U measured from the normal weather 
map is 0.29 radians per day or 17° longitude 
per day. With the value $ = 27 we find that 
m? must be about 15. This corresponds to a 
halfwave-length D of 33 degrees of latitude, 
which agrees well with the value 25 estimated 
from the shape of the continents. 

The solution of the non-homogeneous 
equation (23) is unique when we demand that 
z(x) and its first derivative be continuous 
and periodic with the period 27. By means 
of FOURIER analysis, the solution may be 
written in the form 


z(x) = x42 E h(x)D,(x — x) dx 


where 


n= — > 


Strictly speaking, ©,” has been so chosen 
that the function represented by the integral 
(25) is not the solution to (23), but is rather 
the deviation of the solution from its mean 
valve around the latitude circle. 


Equation (26) may be summed to 


cos s(x — 7 I 
Pol) 25 & ST | ans? OSX 2m; 
Dy (x + 27) = D(x) (27) 


With s = 2.5, this function is represented 
graphically in Fig. 2. 

The formula (25) was applied to a function 
h(x), determined from The Oxford Advanced 
Atlas (J. BARTHOLOMEW 1942) in the following 
way: The height along the latitude circles 
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Fig. 2. The Green’s function ® (x) for the stationary 

perturbation. The dotted curve represents the function 

for the case of no friction, the dash-dotted curve for the 

case of moderate friction, and the dashed curve for the 
case of strong friction. 


40°, 42°, 44°, 46°, 48°, and 50° was estimated 
from the contour lines of the map by taking 
points at 10° intervals of longitude. The 
arithmetic mean of the six values along the 
meridian x was taken to represent h(x). This 
function is represented graphically in Fig. 3. 
With s = 2.5 and x = 0.4 the correspond- 
ing stationary profile of the soo mb surface 
was computed from (25) by approximating 
the integral by a sum in steps of 10° of longi- 
tude. The result is represented by the dotted 
curve in Fig. 4. The agreement with the 
solid curve, representing the normal profile, 
is seen to be poor, and we must again re- 
examine the assumptions. 

We notice that equation (23) gives resonance 
with infinite amplitudes for the solution if s 
is an integer. This is a property of linear 
mechanical systems with no energy dissipa- 
tion. In general one can say that such systems 
are extremely sensitive to small continuously 
acting forces. To remove this artificial sen- 
sitivity it is natural to consider the effect of 
friction. 


The effect of surface friction 


VI. 
With U = 0.29 radians per day, it will take 


an air particle three weeks to travel all around 
the latitude circle. It seems likely that frictional 
damping will be significant after such a long 
time, so that the greater part of the disturbance 
given to the air by a certain mountain will 
probably have disappeared when the air 


returns three weeks later to the same mountain. 
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Fig. 3. Topographical profile at 45° N. The elevations 
represent the mean height of the earth’s surface, averaged 
for the interval 40° N to 50° N, as a function of longitude. 


It is therefore reasonable to assume that we 
have to take friction into account in order to 
explain the stationary perturbations. 

One way of doing this is the following: 
From the theory of the motion in the friction 
layer near the ground, viz., the theory of the 
EKMAN spiral, one can compute the total air 
transport in the friction layer. The formulas for 
the transport are given by BRUNT (1941). The 
total mass transport along the isobars is, in 
the first approximation, proportional to the 
geostrophic transport, and will therefore be 
approximately non-divergent. The mass trans- 
port across the isobars towards lower pressure 
may be written 


H 
D = Fk x 09 Yo 


f zu 


where go, and v,, are the density and geo- 
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Fig. 4. Normal height profile of the 500 mb surface at 

45° N for January together with computed stationary 

profiles for ¢ = o (no friction), for @ = 0.25 (moderate 

friction) and for 0 = 0.50 (strong friction). For purposes 

of comparison the heights are represented as deviations 
from their respective means. 
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strophic wind in the friction layer, k is a unit 
vertical vector, and F is given by 


a (29) 


Here « is the angle between the isobars and 
the surface wind, and K is the “eddy-diffus- 
ivity”. The transport across the isobars means 
that the horizontal pressure forces will do 
work on the air in the friction layer which will 
compensate for the dissipation of kinetic 
energy in this layer by friction. The flow 
across the isobars in the frictional layer will 
produce compensating currents across the 
isobars in the air above the friction layer, 
which in the present case will result in negative 
work being done by the horizontal pressure 
forces and a corresponding decrease in the 
kinetic energy. By this mechanism, the dis- 
sipating effect of surface friction can be transfer- 
red to the air above the friction layer, even 
though the frictional stresses in that part of the 
atmosphere are negligible. 

As an illustrative example, we consider a 
barotropic linear current u along the x-axis 
with constant velocity above the friction layer. 
In the friction layer there is then a mass 


D = + Fou in the y-direction. 
Assuming that the current has a finite lateral 
extent, this means an accumulation of mass 
to the left of the current and a depletion of 
mass to the right, so that the horizontal 
pressure gradient will decrease at all levels. 
In response to this decrease, the air above the 
friction layer will begin to travel across the 
isobars towards higher pressure. This cross- 
isobar current will increase in strength until 
it approximately compensates the mass-trans- 
port towards lower pressure in the friction 
layer; then the pressure field will become 
nearly stationary. This balance is expressed by 


transport 


© 


ll 
f ents =D mn where the in- 
tegral is taken upward from the top, z,, of 


the friction layer. Assuming the cross-isobar 
component v above the friction layer to be 


F 2 4 
constant, we find v = — — y. Since friction 


is supposed to be negligible above the fric- : 


tion layer, the x-component of the equation 


2 ; 
of motion is — = fv. Hence we obtain 


dt 


= —= — Fu (30) 
This equation implies an exponential decrease 
with time of the current velocity, the loga- 
rithmic decrement being F. Adopting the 
values used by BRUNT, viz., « = 22.5°, K = 10 
m?sec-1, and putting f = 10° *sec ?, H = 8,000 
Mm, we obtam F=20 x 10° sec — 017 
day. This means that the velocity decreases 
to r/e times its original value in six days. It 
should be noted however that the decrease will 
be slower in a baroclinic current, where the 
current velocity increases with height. In this 


du 
dt 
where the bar means mean value with respect 


to pressure, and u, is the velocity at the top 
of the friction layer. For instance, assuming 


du du 
Wage 


— — 2—", one obtains a logarithmic decre- 
dt dt 2 


case, we obtain, instead of (30), — = — Fu, 


ment of 0.085 day. 

Similar effects of friction are dealt with by 
A. EINSTEIN in a suggestive essay in his book, 
Mein Weltbild. He points out that friction is 
responsible for the “meridional” circulation 
observed in a cup of tea when the contents are 
set in rotation. He also shows that in a curved 
river, the same frictional effect will cause a 
lateral circulation which is important for the 
explanation of the meandering of the river. 

Since there cannot be any significant accu- 
mulation of mass in the friction layer, we have, 
with good approximation, div D + o,w, = 0, 
where Q and w, are the density and vertical 
velocity at the top of the friction layer (1,000 m 
above the ground, say). Hence, from (28) 


? 


ye ——divD = BR 
0 4 


Thus, ground friction will cause ascending 
motion in regions of cyclonic vorticity near 
the ground and descending motion in regions 
of anticyclonic vorticity. For instance, putting 
H = 8,000 m, F= 2.0 + 10-8 sec-1, and €; =f, 
which means an intense cyclonic vortex, we 


(31) 
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obtain wy = 1.6 cm sec-t. Surface friction 
may therefore be an effect of some importance 
for the maintenance of the typical vertical 
circulations in cyclones and anticyclones. 

The vertical circulation caused by friction 
will give rise to a damping of the perturba- 
tions in the free atmosphere in systems with 
little change of phase with height. For instance, 
in a region with relative cyclonic vorticity 
there will be horizontal convergence in the 
friction layer, a corresponding ascending 
motion at the top of the friction layer, and a 
nearly compensating divergence above the 
friction layer. This divergence gives a decrease 
of the cyclonic vorticity. This effect can 
easily be incorporated into the equivalent- 
barotropic model when we assume that this 
model applies to the air above the friction 
layer. For the vertical velocity ıw, in (6) we 
substitute the value (31). In addition to the 
vertical velocity due to friction, however, 
we have the vertical velocity due to topo- 
graphy, so that the equation of the equivalent 
barotropic model becomes 


This equation was derived under the assump- 
tion that the shape of the streamlines does 
not change with height, and that the horizon- 
tal velocity is of the form v(x, y) A(p). This 
assumption implies that we may write the 
velocity and vorticity near the ground as a 
certain fraction x of the velocity and vorticity 
at the equivalent-barotropic level, i.e., 


Ug, = KU, Vg, = AV, Ce 


We ORE 
HF ) (a sh 3) ate 


Hence (32) becomes 


1 There is no corresponding increase in vorticity in 
the friction layer, because there the frictional forces will 
cause a decrease in vorticity which will approximately 
compensate for the increase due to convergence. 


VII. Explanation of the forced stationary 
perturbations 


We now return to the stationary mountain 
perturbations. By using (34) in linearized form, 
we are able to take friction into account in the 
treatment of the motion. In deriving the 
perturbation equation from (34), the basic 
current is supposed to satisfy the corresponding 
homogeneous equation (h =o). It should 
be noticed that this will not be true if the 
basic current is assumed to be of limited lateral 
extent; then the friction term will not vanish, 
and the basic current velocity will decrease 
exponentially with time. Since such a damping 
is not observed, one has to assume that the 
frictional dissipation is compensated by some 
mechanism not implied in (34) (e.g., probably 
by a slow thermally produced meridional 
circulation), so that as a result the westerlies 
are maintained. This reasoning suggests that 
we may neglect friction in (34) when dealing 
with the basic current, and consider only fric- 
tional effects due to the perturbations. In 
any case, this problem is formally avoided 
by assuming the basic current to be constant. 
Doing so, and assuming the perturbation to 
have a sine dependency on the y coordinate 


02z ’ 2 
mer m?z |, we arrive at the following 
oy? 


perturbation equation for stationary flow, 
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2 
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al 
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where, as in section IV, we set A? = 1. The 
following abbreviations are used: 


E 
3a es m, o= a (36) 


The solution is unique when it is demanded 
that z(x) and its first and second derivatives be 
continuous and have the period 27. It may be 
written 


= 8 J h(a) Px (x — x) da (37) 


where the Green’s function ®,(x) is given by 
+o 
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To evaluate this series, we first determine its 
sum, ®° form = o. We obtain 


—° (x—x) 
D (x) = — = ni C0 [s'(x — 7) — 4], 
A) 
OX S200; 
D (x + 2x) =D (x), 
where 


\/ sinh? @ 7) cos? (sx) + cosh? @ x) sin? (s’x) 


4) 


sin 6 = x sinh (£ x) cos (s’zt), 
cos 6 = 57 coh i. x) sin (sx) 


After having determined D (x), D,(x) was 
calculated by taking 7 terms in the series 
for ®,(x) — D,(x), which converges rapidly. 
The calculation was performed for s = 2.5 
and for o = 0.25 and 0.50, corresponding to 
values of K, the eddy diffusivity, in the ob- 
served range. The associated GREEN’s functions 
are graphed in Fig. 2 together with ®,, the 
GREEN’S function without friction (0 = 0). 
We remark that the curves can be interpreted 
as the stationary height profile caused by an 
infinitely narrow mountain range located at 
x =o. With this interpretation we may con- 
clude from the figure that a narrow mountain 
range will set up a sinusoidal wave with a 
sharp ridge at the mountain. The waves will 
be damped downstream in the case of friction 
and undamped in the case of no friction. 
Since a continent of arbitrary profile may 
be regarded as a sum of infinitesimal ‘‘point”’ 
mountains, we may interpret (37) as ex- 
pressing the general stationary disturbance 
as the sum of the infinitesimal disturbances 
produced by these “point” mountains. 
Using the GREEN’S functions for o = 0.25 
and o = 0.50, the stationary profile of the 
soo mb surface corresponding to the wountain 


profile shown in Fig. 3 was computed from 
(37). The value 0.4 was chosen for the reduc- 
tion factor x; this value seems reasonable from 
wind observations. The resulting stationary 
profiles are shown in Fig. 4. A glance at the 
figure shows that the agreement between the 
computed and observed profiles is very much 
better when friction is taken into account than 
when it is not. Indeed, the agreement in the 
former cases, especially for o = 0.50, is .as 
good as can be expected from a linear theöry. 
It may therefore be considered established that 
the stationary perturbations in the westerlies 
are forced perturbations created by the forced 
ascent of the westerly current over the con- 
tinents and modified by friction. 

The value o = 0.25 corresponds to BRUNT’s 
value K — 10 m? sec-1, whereas o = 0.50 
corresponds to K = 40 m? sec-1. It is known, 
however, that the probable error in the deter- 
minations of K is so high, and the geographical 
variation of K so great, that correspondence of 
K with observed values in order of magnitude 
is all that can be demanded from the theory. 


The revised forecast method 


VII. 


Having arrived at an explanation for the 
persistence of the very long-wave components 
in the observed flow in terms of topography 
and friction, we proceed to incorporate these 
factors into the forecast equation. Accordingly 
we replace equation (6) by (32) or (34). 

Again considering small perturbations on a 
zonal current U, (34) can be linearized to 


Oz 0z oh 
rips ee oe oT 
B Ox Ot cu ù (39) 


where, as in the case of the stationary per- 
turbations, we have assumed a sine variation 
of the perturbation with y given by the rela- 
tion (22). 

In the following, the term ‘‘forced perturba- 
tion” is used to denote any solution of (39) and 
the term “free perturbation” is reserved for a 
solution of the corresponding homogeneous 
equation (h = 0). The most general motion 
satisfying (39) can be regarded as the sum of 
a forced stationary perturbation and a free 
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moving perturbation. Having discussed the 
forced perturbation in section VII, we now 
turn to a consideration of the free perturbation. 

If we consider an elementary wave solution 
of the homogeneous equation of the form 
Z(o,,t).= Ada), we find 


n? + m? 


n? + m= + 2 


n2 — 52 


n? + m? + 22 


+ inU 


UNE" à 


where, as before, s is the stationary wave num- 
ber (24). Thus, a consequence of (39) is that 
the. traveling free perturbations are damped, 
with a logarithmic decrement approximately 
xF. Using the values of x employed in the 
treatment of the stationary perturbation, xF is 
found to be of the order 0.1 day. This 
shows that after a few weeks the traveling 
perturbations would be damped out and only 
the stationary perturbations would remain. 

This damping is, of course, not in agreement 
with the observations. The main reason is 
probably that the traveling perturbations are 
formed as unstable waves, which are able 
to grow in amplitude at the expense of the 
energy of the basic current.! 

As we are here concerned only with 24 
hour changes in the motion, the effect of 
friction will in any case be slight and may be 
ignored. The case is different for the permanent 
motions; there the entire history of an air 
parcel in its travel around the earth must be 
considered, and the effect of friction is decisive. 
We will therefore consider the solution of 
(39) with F =o. 

We have earlier considered the general 
solution as the sum of a free perturbation and 
a stationary forced perturbation, which to- 
gether satisfy the observed initial conditions. 
One may also regard the solution as made up 
of a free perturbation, which alone satisfies the 
observed initial conditions, and the forced 
perturbation produced by the topographical 
forcing of an initially undisturbed zonal 
current. The latter point of view is preferable 
in principle for the motion is then treated 


1 The instability may be attributed to the horizontal 
shear of the current (H. L. Kuo, University of Chicago, 
in a paper to be published in the Journal of Meteorology) 
or to its baroclinicity (J. G. CHARNEY, 1947). In both 
cases the very long wave forced perturbations are stable 
and are therefore subject to the full damping effect of 
friction. 


purely as an initial value problem, or as 
RICHARDSON called it, a ‘‘marching problem”. 
This method makes it unnecessary to determine 
the stationary perturbation, which is a problem 
“in the large”, and enables us to look at the 
forecasting problem as a “local” problem, the 
solution of which depends only on what is 
happening within a limited region. 

Adopting the latter point of view we 
obtain the solution of (39) (with F =o) by 
putting 


+ © ++ © 


ae.) = (er, h(x) = > ap, er 


n= — oo n=— @ 


Substituting these expressions into the differ- 
ential equation, we obtain a differential equa- 
tion of the first order for c,(¢). Solving this, 
and adjusting the constants so as to satisfy the 
initial condition, we finally get 


zx + Ut, t) =z (x; 0) + 


+ f z(a,0) Le (x — «x, t) da + 


+ x 22 yp (x) J (x —«a, t) de, 


where 
ace inbt 
Lats, iy = >: (er a :) ed) 
n=—o 
fe ind ar 
I(x, = I D (au cu ) aa (42) 
n = — CO 
and 
a — m? + 2 = 2 52, b = B+ 22U, 


be =F ne 


(43) 


Let us consider the first integral in (40). 
This is the solution of the homogeneous 
equation satisfying the initial condition. The 
influence function Lis (x, f) is essentially the 
same as dealt with in section IV, the only 
difference being that A? is replaced by a?, since 


so 


Io? (x,1) 


— 
-180 -I50 -I20 


Fig. 5. The influence function Ia? (x, 1) for = 2.5 
and a?= 10. This function represents the 24 hour 
height change of the soo mb surface, relative to a co- 
ordinate system x traveling with the mean current, pro- 
duced by a concentrated initial disturbance at x = 0. 


we now consider perturbations with finite 
lateral extent. In the study of the stationary 
perturbations we found that m? should be about 
15; consequently, we should expect a? = 18 
to be a proper value. The influence function 
was computed in the same way as shown in 
section IV, for a time interval of one day 
(j= Band fon 6,30, 14,and. 18.) Lhe 
results are given in table I, and, in the case of 
a? = 10, also in Fig. 5. There turned out to 
be very little difference in the forecasts made 
with the influence function for a? = 10, a? = 
= 14, and a? = 18. The value a? = 10 used in 
the following was based on a previous estimate 
of m which differs somewhat from the value 
given in section V. Because of the relative 
insensitivity of the forecast to changes in m it 
was not found necessary to revise the results. 

It will be seen that the influence function 
for one day is almost zero in the greater part 
of the region (— 180° to 180°), and has an 
appreciable value only in the vicinity of the 
origin. It is therefore not necessary in (40) to 
integrate all around the latitude circle; it is 
sufficient to take the integral over a certain 
influence region around the point x. This 
procedure is in agreement with group velocity 
considerations (CHARNEY, NWP) and means 
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that it is not necessary to know the data for 
the entire world in order to make a 24 hour 
forecast in a small region. In this respect, it 
seems preferable to use the influence function 
with a? = 18 since then the influence region 
would be smaller. 

The last term in the forecast formula (40) is 
a solution of the non-homogeneous equation 
and represents the disturbance produced: by 
the continents after the time f in an initially 
undisturbed current: This part is therefore 
independent of the initial profile z (x, 0) and 
can be determined once and for all as long as 
s and U remain the same. 

In order to compute the GrEEN’s function 
J(x, #1, we write (42) in the form 


+ 00 + 00 ; 
T ein(x+Ut) ein 
Nn) Eee an ms eee 
Je 27 ; ‚n?—s? Haze 


n = — © n=— CoO 


The first two series can be summed exactly, 
and the last series converges so fast that sum- 
mation from — 20 to + 20 gives sufficient 
accuracy. In this way, the function J(x, 1) was 
computed for a? = 10 and U =0.3s radians 
per day or 20° longitude per day. The result 
is shown in Fig. 6. The shape of the, curve 
shows that there is a finite influence region 
for the mountain influence also, which means 
that the forecast is not influenced by mountains 
very far away. The width of the influence 
region is about the same as for the influence 
function I,, (x, 1) for free perturbations. The 
function J (x, 1) may be interpreted as the 
disturbance produced in an initially undisturbed 
current after 24 hours by an infinitely narrow 
mountain range situated at x =— 20°. The 
reason for the two discontinuities in the curve 
is that the air between the discontinuities has 
passed over the mountain during the period, 


! The method given for calculating J(x, f) fails when 
s is an integer, corresponding to resonance for the sta- 
tionary perturbations. However, J(x, ft) is finite in this 
case too and can be determined by a different method. 
That J(x, f) must be finite is obvious from energy con- 
siderations. 
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Fig. 6. The Green’s function J(x,1) for the 24 hour 

height change of the soo mb surface at latitude 45° pro- 

duced by the forced motion of the current over an 
irregular surface. 


whereas the air outside the discontinuities has 
not. 

The corresponding 24 hour disturbance 
produced by the actual mountain profile is 
represented by the last term of (40). This 
term was computed with x = 0.4 for the 
mountain profile h(x) shown in Fig. 3. The 
result is shown in Fig. 7. 

It is interesting to note that the Rocky 
Mountains act approximately as an infinitely 
narrow mountain range. The computation 
shows that the topographical influence on a 
24 hour change in the pressure field is small. 
In fact, in applying (40) to actual weather 
situations, it was found not worth while to 
take the topographical term into account for 
a 24 hour forecast, even though it would have 
improved the forecast slightly. It should be 
remembered, however, that the effect of 
topography will be larger for a longer forecast 

eriod. 

Without the effect of topography, and 
with a limited influence region, the prognostic 
formula (40) becomes 


Ar HUN (60). + 


X + X 


+ VE TC? 0) La? 


x—ı 


(w—a, 1) de. (44) 


This formula was tried out on some actual 
situations taken from the U.S. Air Weather 
Service map series Historical Weather Maps for 
January, 1946. The function z (x, f) was identi- 
fied with the height of the soo mb surface at 
45° latitude. The zonal wind U was taken to 
be 20° per day, and the value 4? = 10 was 
used in the influence function. The integral was 
taken over an influence region of 120° longi- 
tude (x, = 30°, x, = 90°) to insure accuracy, 
although in practice it would be possible to 


reduce the total interval to 90° or less. It was 
evaluated by means of Srimpson’s formula 
using steps of 10° of longitude. The formula 
actually used was 


220% 2) = 2% 0) 
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+ 2 Ay z(x + N-10°, 0). 


N=—3 


(45) 


The values of the coefficients Ay are given 
in table II for a? = 10, 14 and 18. In this 
way, a forecast for a region of 90°, say, can 
easily be computed in a half hour. Since the 
influence function has the mean value zero, 
we may subtract any constant from the height 
z in computing the sum. In order to prevent 
having too large a value for z, this quantity 
was reckoned from the level 18,000 feet. 


IX. Discussion of results 


The results of 6 forecasts are presented in 
Figs. 8—13. The upper curve in each diagram 
shows the initial profile of the soo mb surface 
along the 45° latitude circle. This curve was 
obtained from the soo mb maps by reading 
off the heights at every 10° of longitude to 
within 10 feet. The two curves in the lower 
part of each diagram show the computed 
(dashed curve) and observed (solid curve) 
local changes in the following 24 hours. 

In judging the results, it should be borne in 
mind that, owing to the scarcity of aerological 
data, the soo mb maps involve errors which 
often may be of the order of 200 feet. Since 
the data on which the analysis is based are most 
extensive in the region 130° W—80° E, the ana- 
lysis will presumably be most reliable in this 
region. This is also the region where the 
method gave the best results on the average. 


-5 ! 
180° 120° 60° 00° 60° 


120° 180' 

W E 

Fig. 7. The 24 hour height change of the soo mb 

surface at 45° N, produced by the motion of an initially 
straight zonal current over the continents. 
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Fig. 10. likes agi 
Figs. 8—13. In each figure, the upper curve represents the observed height profile of the soo mb surface at 


45° N for the date indicated in the upper right hand corner. The lower continuous curve gives the observed 
change in the following 24 hours, and the dashed curve the predicted change. 
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The diagram for January 13 shows an obvious 
example of a disagreement between the com- 
puted and observed changes due to erroneous 
analysis. The map for January 14 shows a 
tremendous filling of the entire Pacific trough, 
which is probably fictitious since there is 
very little data to substantiate it; moreover, 
it deepened again the next day. The large 
“observed” pressure rise in the region 90° E 
to 180° E, which is shown in the diagram for 
January 13, is therefore entirely wrong; in 
this case, the computed forecast was probably 
much better than the analysis for January 14. 

Besides the errors due to erroneous analysis, 
we have to consider errors due to possible 
non-fulfillment of the assumptions underlying 
the forecast. The forecast equation was derived 
for small perturbations on a zonal current, and 
one cannot expect good results when the 
amplitude grows too large. Furthermore, the 
perturbations were assumed to vary in a 
certain regular manner in the y-direction; 
this assumption made it possible to represent 
the two-dimensional pressure field by a profile 
along one latitude circle and thereby to treat 
the problem as a one-dimensional one. When 
the actual perturbations do not show such a 
regular variation in the y-direction, there is 
always the possibility of an advection of 
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vorticity from the north or south which 
cannot be predicted from the shape of the 
profile at one latitude. The strong deepening 
that took place in the Atlantic (40° W) from 
January 9 to 10, and which the method failed 
to predict, may be due, at least partly, to the 
advection of strong cyclonic vorticity from the 
north. 

Errors may also result from an improper 
choice of the current velocity U. Even though 
the derivation of the forecasting formula was 
made for a constant U along the latitude 
circle, it is possible to a certain extent to 
account for slow variations of U with longi- 
tude. This is so because the forecast is independ- 
ent of the zonal current outside the influence 
region. The best way to choose U has to be 
determined more or less empirically. The 
authors’ choice of a constant U equal to 20° per 
day for all six cases is of course a very rough 
estimate and could probably be improved. 

The errors mentioned so far are due to the 
assumption that the motion can be dealt 
with in one dimension, using the linearized 
equation. These errors will therefore be eli- 
minated when it is possible to integrate the 
non-linear equation for two horizontal dimen- 
sions with arbitrary initial conditions. (This 
will probably be possible in the near future, 
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but the computation will be very elaborate 
and will require high-speed and high-capacity 
computing machines.) It seems very likely 
that a computation by means of the non- 
linear, two-dimensional equivalent-barotropic 
equation would improve the forecast consider- 
ably. Since some success is obtained by using 
the simplified one-dimensional linearized equa- 
tion, we may be justified in concluding that 
the quasi-geostrophic approximation and the 
equivalent barotropic model are legitimate 
approximations that render sufficient accuracy 
to make them useful tools in numerical 
forecasting. 


X. General remarks on the stationary 
perturbations 


It has for some time been recognized that the 
quasi-stationary perturbations of the atmosphere 
are caused by geographically fixed perturbing 
forces, but the exact nature of these forces has 
not been well understood. RossBy (1939) was 
able to show that their effect would be to 
excite a system of free perturbations, whose 
dimensions are in qualitative accord with 
the observed large-scale stationary perturba- 
tions in middle latitudes. The establishment 
of the continents as the primary perturbing 
influence enables us to fill the blanks in the 
picture and to determine not only the mean 
wave lengths of the system but also the posi- 
tion of the individual troughs and wedges. 
These are not anchored by solenoid fields 
along coast lines, as has been sometimes 
thought, but are determined by a world-wide 
adjustment of the zonal flow to the shape of 
the continents. 


~ 


The relationship between zonal index and 
mean circulation gives the possibility of 
explaining the general longitudinal variation 
of climate as a function of the zonal index 
for many surface climatic factors are cor- 
related to the upper flow, which is in turn 
determined by the index. Thus, because 
of the horizontal barotropy of the large- 
scale stationary perturbations, we may say 
that the mean surface temperature is correlated 
to the mean soo mb height — except in those 


cases where the surface temperature is not 


representative of the free air temperature. 
Similarly, since the mean position of the polar 
front parallels the upper current, cyclone 
tracks are also determined by a knowledge of 
the stationary circulation pattern. 

In this way, we have established a physical 
relationship between the zonal index and the 
main weather type, a relationship which has 
hitherto been considered mainly from a sta- 
tistical point of view. This result has far reach- 
ing consequences. If it should prove possible to 
forecast the variations of the zonal index, it 
would also be possible to make at least an 
elementary forecast of the change in world 
weather. The consequences of the relation- 
ship for longrange weather forecasting, as 
well as for the study of climatic changes, are 
obvious. 
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Recordinds of Cosmic Radiation 


By K. G. MALMFORS 
The Royal Institute of Technology, Stockholm 


(Manuscript received 22 April 1949). 


Abstract: 


The main purpose of the registrations reported here has been to study the diurnal 
variations of C. R. and especially to investigate how the amplitude and the phase of the 
variation depend on the direction of incidence. Since statistical Auctuations will dominate 
over the significant ones unless a sufficiently large number of coincidences is counted, it 
has been necessary to use an apparatus with a high counting speed. C. R. has been 
recorded continuously in Stockholm (magn. lat. 58°) since 1. 10. 1947 by means of large 
counters in coincidence circuit. Recordings have been made in three different directions 
simultaneously, and the total number of particles counted amounts to 10°. The material 
has been analyzed with respect to diurnal variation. This analysis definitely proves that 
the variation is not the same in different directions. A term called C. R. activity has 
been introduced for the purpose of classifying the registrations according to the magni- 
tude of the variation for each day. An analysis shows that large activity is very often 

observed during consecutive days. 


Introduction 


The intensity of C.R. recorded at the 
earth’s surface undergoes only small varia- 
tions. The largest ones are due to atmospheric 
influences. Since the absorption of C. R. in 
the atmosphere depends on the amount of 
matter traversed the intensity of the radiation 
is closely correlated to changes of barometric 
pressure. The coefficient can be accurately 
determined and consequently it is easy to make 
a correction for such variations. On the other 
hand the intensity shows some correlation 
with the air mass distribution as was first 
pointed out by BLACKETT (1938), who inter- 

reted this temperature effect as being due 
to the fact that, when the atmospheric tem- 
perature is high, the mesons originate from a 
higher level, which leads to a greater pro- 
bability for their decaying before reaching 
the earth’s surface. The effect has recently 
been studied by DUPERIER (1945, 1948), who 
correlated the intensity of C. R. with the 
height of the too-millibar surface determined 


from balloon sondations. A correlation cocffi- 
cient as high as 0.67 was obtained and the de- 
pendance was determined at 5.4 %/km. 

Cosmic ray intensity is very intimately 
associated with magnetic storms and, even 
though this connexion shows many complex 
features, it seems probable that investigations 
following this line will give essential know- 
ledge of the nature of the magnetic storms. 
The analysis of the present registrations in 
this respect will be reserved for a later publica- 
tion. 

The diurnal variations of C. R. with respect 
to solar and siderial time have been studied 
by many investigators. Recent results con- 
cerning variation with solar time show good 
agreement and indicate that the intensity is 
some few tenths of a percent higher in the 
afternoon than in the morning. On the other 
hand, the existence of some variation with 
siderial time has not been established with 
any certainty. Results indicating such a varia- 


56 


tion may probably be interpreted as being 
due to seasonal changes of the variation with 
solar time. 


Apparatus 


The apparatus described here is mainly an 
enlarged copy of an apparatus used for con- 
tinuous registrations during the period 1939— 
1943 (ALFVÉN and MALMFORS 1943). 

The counters used are made of brass tubes 
(diam. 15 cm, length so cm) closed at the 
ends by plane brass plates. The axial wire 
(Wo, 8 =0.2 mm) is suspended by using 
glass metal scals. Filled with carefully dried 
argon-alcohol mixture (10 mm C,H,OH, 
35 mm Ar) the counters operate satisfactorily 
at about I 400 volts. The length of the plateau 
is usually about 200 volts. The high voltage 
for the tubes is supplied from dry batteries in 
order to ensure the highest degree of con- 
stancy. 

Ten tubes arranged as shown in Fig. ı 
make it possible to record the intensity in 
twelve coincidence sets simultaneously, (4 
channels recording the intensity from the 
zenith, 4 channels that from a cone inclined 
to the north and 4 channels the intensity 
from the south). The number of coincidences 
in each of these channels is recorded separately 
by an ordinary coincidence amplifier with 
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Rossi circuit, scale-of-four and a mechanical 
recorder fed by a thyratron. 

The anode- and grid voltages and the fila- 
ment currents of the amplifiers are stabilized 
in such a manner that no appreciable varia- 
tions can arise from fluctuations of the main. 

By counting the intensity from each direc- 
tion in four independent channels valuable 
control is obtained. If one of the counters or 
one of the amplifiers does not operate satis- 


factorily this is easily detected by comparing 


the four channels. In such a case three of them 
will show good parallelism while the fourth 
varies irregularly and can be repaired without 
interrupting the counting in the three un- 
disturbed channels. 

The recorders are photographed once every 
hour by means of a signal clock and a relay 
system. 

The apparatus is placed in a little cottage 
a hundred meters away from other buildings. 
The temperature in the room is kept constant 
to within some tenths of a degree. On days 
when temperature variations exceeded o.5° C 
the registrations have not been subjected to 
analysis. 


Representation of data 


Tab. I, which is a reproduction of the photo- 
graphic recordings, shows how the counting 
capacity of the apparatus 115 000 coincidences/ 
hour, is distributed between the different chan- 
nels. (It should be noted that, owing to the 
scale-of-four, the real number of coincidences 
is four times as great as that derived from the 


28cm photographs.) The results given here com- 

prise registrations during the time 1/10 1947— 

8/1 1949, the total number of coincidences 

amounting to about 10%. The first step in 

< BER > the treatment of the material has been to 

plot the 24-hour values of each channel, in 

Fig. 1. Arrangement of the counters. order to control the counters and the ampli- 

Table I 

RQ es 

3/7 1948 Z N S 
ore? 8774 | 6662 | 6235 | 5151 | 4000 | 1472 | 0602 6556 | 8793 | 7893 | 00 

2 S 34 | 6933 | 6213 

022° 8775 | 9636 | 9296 | 8166 | 7061 | 3382 | 2447 | 8508 7013.1.97207 | 7831.21 8773 |. onsu 

030° 8776 | 2581 | 2458 | 1223 | oror | 5289 4375 | 0465 | 3201 | 1567 | 3635 | 0604 | 0126 

04 8777 | 5527 | 5528 | 4296 | 3128 | 7189 | 6272 | 2395 | 5351 | 3410 | 5458 | 2434 | 2099 
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Table II 
3/7 1948 Z | N | S Bar. 7 | N | S Z | N | S | 
| 
022 24167 15797 15010 758.6 24089 15746 14962 2.5 m7 3.2 
200 24315 15907 14859 758.6 24237 15850 14811 3.1 2.4 2.2 
490 24242 15863 14904 758.6 24164 15812 14856 2.8 D 


fiers. If one of the channels has happened to 
vary irregularly, it has been excluded. The 
remaining channels have been summed up and 
2-hour values have been calculated. The 


/947 December 1949 Januory 
BR LE EN EI CERT OTE TE CR CE CE 


first three columns of Tab. II contain values 
obtained in this way. Barometric data have 
been kindly supplied by The Meteorological 
Station, Stockholm, and during the whole 
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Fic. 2. Records of C. R. intensity. The upper, middle and lower curves correspond to Z, N and S directions 
a respectively. The distance between the lines is 2%. 
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Table III 
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Disturbed days 


Quiet days 
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series of registrations the barometric coeflicient 
— 0.23 %/mm Hg has been used. No cor- 
rection has been applied with regard to the 
atmospheric temperature effect. After reduc- 
tion to standard pressure (760 mm Hg) the 
material takes the form shown in the fol- 
lowing three columns in Tab. II. Owing to 
the fact that one or more of the channels 
have sometimes had to be excluded, the 
corrected 2-hour values have been represented 
(in the last three columns) by the percentage 
deviation from a mean value. From the points 
of change which have always been made to 
occur at midnight, new mean values have 
been calculated. This has been done in such a 
way that the mean values of the curves for 
the preceding and the following day change 
in the same proportion as the mean values 
of the continuating channels. 

An analysis of the material with respect to 
atmospheric influences will be published by 
Dr. LINDHOLM. 

In Fig. 2 some parts of the registration series 
have been reproduced. The upper, middle 
and lower curves correspond respectively to 
the Z, N, and S directions. The distance 
between the horizontal lines is 2%. The 
period 10/3—16/3 1948 is one of large magnetic 
activity including a severe magnetic storm 


! The author will be pleased to send copies of the 
whole registration series to those interested. 
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beginning on March 15. Since this represen- 
tation of the data has been chosen for the 
analysis of short time variations, such channels 
have also been included which slowly and 
regularly decrease with time, which some- 
times has been the case. No great importance 
should therefore be attached to apparent 
variations ranging over periods as long as 
one month or more. The long time variations 
are more accurately obtained by selecting 
only such channels as give perfectly parallel 


curves. 


The diurnal variation 


A systematic diurnal variation can only 
occasionally be distinguished from the sta- 
tistical fluctuations when studying the curves 
of Fig. 2. Significant results can only be 
obtained by calculating a mean variation 
covering several days. Each point in the 
harmonic dial diagrams of Fig. 3 has been 
obtained by summing up the results of six 
consecutive days. The mean values are re- 
presented by the square points. Fig. 4 shows 
curves obtained from registrations during in- 
tervals of 100 consecutive days. The similarity 
of the curves confirms the reproducibility of 
the variation. From any representation of the 
data as given above it is obvious that the 
variation is not the same when measured in 
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Fig. 3. Harmonic dial diagrams of diurnal variation. Each point corresponds to a period of 6 days. The 
square points represent mean values. 


different directions. The curves for the difference, 
N minus S, have about the same amplitude 
as the individual curves. Similar results have 
been obtained from earlier registrations by 
KOLHÔRSTER (1941) and by ALFVÉN and 
Matmrors (1943). 

The harmonic diagrams of Fig. 3 show 
that the dispersion of the points for N—S is less 
than that for the original directions, a fact which 
may indicate that the phase difference between 
the curves for N and S is due to different 
deflection angles in the earth’s magnetic field. 

An analysis based upon the harmonic 
coefficients, calculated for each day, shows 
that the amplitude of the variation is not 
distributed at random. On several occasions 
large amplitues occur during consecutive days, 
as is also illustrated in Fig. 2. During the period 
24—27/12 1947 a distinct variation, lasting for 


some days, is measured in all three directions. 

Since this effect will probably be of im- 
portance for the interpretation of the varia- 
tions, the term “Activity of Cosmic Radia- 
tion» has been introduced, referring to the 
magnitude of the variation during individual 
days. Analogous with the procedure adopted 
for the treatment of magnetic observations, 
the material has been classified into quiet and 
disturbed days. In the case of the present re- 
gistrations this has been effected in the fol- 
lowing way. To the amplitude intervals 
O—0.25, 0.25—0.40, 0.40—0.60, >0.60 % have 
been assigned the indices 0, 1, 2, 3 respectively. 
In this way the C. R. activity for each day is 
given by a total index which is the sum of 
those for the Z, N and S directions. For 72 
days which are called disturbed days this 
total index is 6 or more. 107 days have the 
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Fig. 4. Curves showing the diurnal variation during periods of about 100 days. 
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Fig. 5. Harmonic dial diagram, showing mean values 
of diurnal variation during days grouped together 
according to C. R. activity. 


index $ or 4. 133 days have the index 3 or 2 
and, finally, 87 quiet days have the index 
0 or T. 

Tab. IV which gives the number of cases 
where disturbed and quiet days occur conse- 
cutively makes it evident that the amplitu- 
des, as mentioned above, are not distributed 
at random. 

Fig. 5 shows the mean values of the diurnal 
variation for the four groups of days. This 
analysis gives the result that the phase of the 
variation is nearly the same during quiet and 
disturbed days, and indicates that even for 
quiet days there is a phase difference between 
the different directions. 


Table IV 
Consecutive days 
Total 
oP ee 
Disturbed ..| 29 12 5 I 72 
Outre: 36 13 3 4 87 
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As yet there is no generally accepted inter- 
pretation of the diurnal variation of C. R. 
From the fact that the variation is not the same 
in different directions it is obvious that atmo- 
spheric effects alone cannot be responsible 
for the variation. Furthermore, the heating 
of the atmosphere during the afternoon seems 
likely to give a variation which reaches its 
maximum during the early hours of the mor- 
ning. It seems likely that the observed diurnal 
variation must be interpreted as being due to 
some radiation anisotropy which is closely 
related to the sun. The fact that the amplitude 
of the variation is very often large during 
consecutive days indicates that the strength of 
this anisotropy varies from time to time. 

As a consequence of new ideas on the gene- 


ration of C. R. (ALFVÉN, RICHTMEYER, TELLER, 
1949) it seems likely that the variation may be 
associated with the electrodynamical conditions 
prevailing in the region around the sun. 
According to ALFVÉN, the particle beams 
emitted from the sun during magnetic storms 
are associated with electric fields with potential 
differences as high as 10° volts. ALFVÉN has 
interpreted the change of C. R. intensity 
during a magnetic storm as being due to such 
electric fields (1946). If this interpretation is 
correct the same argument would apply to the 
diurnal variation, and it ought to be possible 
to find some relation between the diurnal 
variation of C. R.and the structure of the electric 
fields as revealed by magnetic disturbances 
and world-wide changes of C. R. intensity. 
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The Trough-and-Ridge diagram 


By ERNEST HOVMOLLER 


Swedish Meteorological and Hydrological Institute, Stockholm 


Abstract. 


A discussion is given concerning the construction and interpretation of a diagram which 
shows longitude and intensity of troughs and ridges at the soo mb level in middle 
latitudes as a function of time. 


1. Troughs and ridges as characteristic 
features of upper air charts 


As. synoptic-aerological observation material, 
during the 1930's, became sufficiently copious to 
permit the regular construction of upper air charts, 
it was soon found that these charts were not of such 
complex character as the surface maps often are. 
Many of the minor singularities of surface maps 
(and surface weather) are due to local orographic 
effects, and it is evident that such irregularities can 
not extend to any great height. Besides, several 
(though not all) minor irregularities of non-oro- 
graphic character were not found to exist — in 
some cases one might say: were found not to exist 
— in the upper air. 

Another characteristic feature of the upper air 
charts is the preponderance of a northward slope of 
the isobaric surfaces in middle latitudes. The pre- 
dominance of the westerly winds which is observed 
near the surface in both hemispheres between la- 
titudes 35 and 65°, is still more pronounced in the 
middle and upper troposphere and the lower stra- 
tosphere. The great cyclonic and anticyclonic sys- 
tems of the surface layers are not smoothed out 
on the soo mb charts, with the exception of some 
anticyclones in which the air is extremely cold; but 
in many cases, the cyclones appear as troughs open 
to the north and the anticyclones as ridges open to 
the south. 

The complex nature of surface charts is a great 
hindrance to any systematical classification of these 
charts and to any other effort towards a rational 
prognostic method. As far as upper air charts — and 
more specially the soo mb charts — are concerned, 
the significant variable factors are not quite sonumer- 


ous: the location and the velocity of the strongest 
westerly flow; the longitude, orientation, and inten- 
sity of troughs and ridges; the location, shape, and 
intensity of cut-off vortices. In the present paper, we 
are mainly concerned with the troughs and ridges: 
their position (in terms of longitude), their intensity 
and intensity variations, and their displacements. 
For our purpose, it is, within certain limits, justified 
to treat the cut-off vortices as over-developments, 
so to say, of troughs and ridges. 

The troughs appearing on 500 mb charts are of 
many different types. Some of them are the insigni- 
ficant attendants of moderately developed warm- 
sector cyclones. Others are minor troughs existing 
in the greater part of the troposphere but of slight 
intensity in all levels. Troughs of moderate intensity 
are often accompanied by well-developed wave 
cyclones in the first phase of occlusion, or of a cy- 
clone which is already filling; others are associated 
with troughs in the lower troposphere and: if so, 
they are in most cases of a non-frontal character 
(so called ‘polar troughs”). The major troughs, 
normally existing to anumber of 2—6 in the northern 
hemisphere, are always cold troughs; very often the 
amplitude of the isotherms is even greater than that 
of the contour lines. 

Minor and moderate troughs move, with only 
few exceptions, eastwards; in such cases where the 
upper trough is associated with a cyclone or a trough 
near the earth’s surface the velocity of the upper 
system is approximately the same as that of the 
system at the surface, though minor differences 
sometimes may occur. The propagation of. polar 
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troughs is generally rather slow; evidently, the air 
moves into such troughs from the west and leaves 
them again towards the east. 

The major troughs of the upper air are nearly 
always slow-moving as compared with wave cy- 
clones; in some cases, they are even stationary or 
have a sligth retrograde movement. 

The life-history of an individual trough is not 
always simple and clear cut. The ideal case, as it 
might be called, of the growth of a minor trough 
into a major one and the subsequent continuous 
decline and disappearance, is not often seen; in most 
cases, We see on our charts rather complicated pro- 
cesses of amalgamation and partition. It should be 
pointed out, however, that many major troughs and 
some troughs of moderate development may be 
followed without difficulty throughout several 
days or even weeks. 

As for the ridges, we find an almost equally great 
number of variations in type and life-history. Ridges 
of maximum development are in most cases almost 
stationary or even slightly retrograde; normally they 
exist for a week or more, thus exerting a ‘‘blocking 
effect” upon wave trains approaching from the west. 


2. Construction of trough-and-ridge diagrams 


The existence of a trough at a certain longitude 
within the belt of latitude which we are considering, 
comes out very clearly if we plot the mean height of 
the soo mb level within that latitudinal belt as a 
function of longitude. Major and minor troughs, 
troughs with steep or weak “‘slopes”’ to the east and 
west, symmetrical and asymmetrical troughs are 
easily discriminated by this graphical representation. 

If we want to follow the movement of a certain 
trough or ridge during several days, a different 
method of representing the observed facts must be 
used. We may simply plot the longitude of the 
trough’ (as found by the method described above or 
by any analogous techniques) as a function of time. 

These two simple ways of illustrating charac- 
teristic features of troughs and ridges may be com- 
bined in a graph showing the mean height of e. g. the 
soo mb level between fixed latitudes as a function 
of longitude and time. Such diagrams have been 
constructed for the period October, 1945, to April, 
1946. 

The material on which the diagrams were based is 
the series of Historical Weather Maps published by 
the Headquarters of the U. S. (Headquarters, Air 
Weather Service: Northern Hemisphere Historical 
Weather Maps). From these maps the height of the 


500 mb surface was taken at the intersection points 
between every fifth degree of longitude and the 
parallels of 35, 40, 45, 50, 55, and 60° N. The choice 
of parallels is to some extent arbitrary. Although the 
annual geographical variation of the belt of maxi- 
mum zonal wind might be taken as a reason to 
choose different latitudes for different months, it 
was considered more advantageous to use the same 
belt for the whole period covered by the investiga- 
tion. The lower limit, 35°, is in the northern border 
region of the subtropical belt with high pressure and 
a rather small slope of the isobaric surfaces; the upper 
limit, 60°, is generally to the north of the belt in which 
the zonal wind in the soo mb level has its maximum 
intensity. Besides that, the aerological material 
north of 60° during the period concerned was too 
scarce to secure a sufficiently correct analysis; in 
fact, the analysis over the greater part of Asia and 
the Pacific ocean is rather tentative even for the 
belt between 35 and 60° N. 

Taking the average of the six values obtained from 
each map to represent every sth degree of longitude 
between latitudes 35 and 60° N (or, as one might 
rather say, between 321/: and 62:/2° N), we get a 
sufficiently good approximation to the actual mean 
value for this part of the meridional circle, provided, 
of course, that the maps themselves are satisfactory. 
By computing such average values for every fifth 
degree of longitude, we are sure to get a rather detailed 
picture of the profile along the zonal belt. The time 
interval between two consecutive maps, 24 hours, 
is, however, rather too long. With an interval of 
12 hours instead of that, the identification of minor 
troughs and ridges should probably be less am- 
biguous. 

As far as the regions with relatively ample obser- 
vation material are concerned, there is no reason to 
smooth the average values plotted in the diagram, 
neither before plotting them nor in drawing the 
isolines. It was thought that little should be gained 
and some unwished subjectivity introduced by 
smoothing the more doubtful values between appr. 
60° Eand 150° W, hence, for the analysis, all values 
were taken to be correct. That means, of course, 
that as far as this longitude interval is concerned 
little or no confidence could be taken in minor 
irregularities occurring in the analyzed diagrams. 

The ideal way of representing the material would 
be to fold the diagram around a cylindre, the time 
scale being parallel to the axis. To reduce the disad- 
vantages of a discontinuity along a certain meridional 
circle, the strip between 120 and 180° E has been 
placed left and right in the diagram. 
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3. Analysis of the trough-and-ridge diagram 
for Nov. 1945. 


The diagram for Nov. 1945 (Fig. 1) shows at a 
glance that an almost stationary major ridge prevailed 
throughout the month near longitude 0°. To the 
west of this blocking ridge, several troughs and rid- 
ges are moving east during the month. Some of the 
troughs seem to penetrate into the stationary ridge 
but this is always rebuilt very soon. The phase 
velocities of troughs and ridges are given by the 
orientations of the “bottom line” of the troughs 
and the ‘‘peak line” of the ridges. They can be 
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Pig. 1. Geopotential (given in dynamical decametres) 
of the soo mb level as a function of longitude and 
time. The values on which the diagram is based are 
average values of geopotentials at 6 different points on 
the same meridian (at latitudes 35, 40, 45, 50, 55 and 
60° N respectively). Areas of high geopotential values 
(i. e. ridges) are shown by horizontal hatching, areas of 
low values (troughs) by vertical hatching. The slanting 
straight lines indicate a succession of maximum develop- 
ment of troughs and ridges which is of the general 
character required by the dispersion theory of atmos- 
pheric waves (see text). 


measured directly in terms of degrees longitude per 
24 hours; for the interval between 180° W and 
the stationary ridge they vary between 5 and 15° 
long. per day. No systematical difference is found 
between the velocities of major and minor troughs 
in this part of the diagram. To the east of the station- 
ary ridge, conditions seem to be complex and ill- 
defined but the main feature may be said to be an 
almost stationary trough, situated at appr. 140° E 
during the greater part of the month. 

It is by no means easy, even in this diagram where 
zonal inequalities have been eliminated, to count 
the number of troughs and ridges at any given time. 
But if we look upon the month as a whole, weiget 
the impression that only two major troughs are 
present, one near the east coast of Asia, another 
(which is more variable in position and intensity) 
near the east coast of North America. The ridges are 
one over western Europe or the eastern Atlantic, 
and one (fluctuating rather much) over the Pacific. 

Looking more closely upon the conditions of the 
western hemisphere as shown by the diagram, we 
will find a remarkable rhytm in the intensity of 
troughs and ridges: An intensity maximum of a 
trough is in several cases followed after one to three 
days by a maximum of the ridge nearest ahead of 
(to the east of) the trough. Good examples of this 
phenomenon are found, above all, during the first 
ten days of the month. Although there are some 
cases where the rule quoted above does not hold 
good, it would seem that the phenomenon is a real 
one and not a matter of chance. In fact, it is essenti- 
ally in agreement with the theory of dispersion 
introduced in meteorology by Rosssy (1945). 
According to this theory, we should have (at 45° N, 
which is nearly in the midst of the latitudinal belt 
with which we are concerned) 
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where c is the phase velocity, cg the group velocity, 
U the velocity of the basic current, and L the 
wave length. (The unit of length used here is the 
length of 1 degree of longitude at the latitude con- 
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cerned, i. e. 1IT--— km, and the time unit is 24 
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hours.) In these formulas, the value of U is most 
difficult to fix. By inserting empirical values (taken 
from the figures) 

c = 10° long./day, cg = 28° long./day, L = 65° 
long., we get from the first equation 


U = 22° long./day = 70 km/t 
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and from the second 
U = 16° long./day ~ so km/t. 


These values of U are quite reasonable and not too 
different from one another. It should be remembered 
that the elementary barotropic theory is based on the 
assumption that the perturbations (the additional 
velocities) are very small as compared with the velo- 
city of the basic zonal current. This assumption, 
which means a serious limitation to the application 
of the theory, is in the present case definitely not 
in good agreement with the observed facts. One 
might, however, conclude that the dispersion pheno- 
menon resulting in a group velocity different from 
(greater than) the phase velocity is not confined to 
such cases where the perturbation amplitudes are 
small. 


4. Brief description of diagrams for other 
months 
(Oct. 1945, Dec. 1945 —April 1946). 


Space does not permit us to publish the other dia- 
grams of the type of fig. ı which have been con- 
structed but we shall give a short description of 
their main features and a figure showing the “wave 

LE = 
profile” of each month as a whole (Fig. 2). 


Fig. 2. Monthly mean profiles of the soo mb level in 
middle latitudes (32,5—62,5° N). 


The number of major troughs is generally 2 or 3. 
The most permanent feature of the whole period 
was a trough near 140° E, i. e. in the vicinity of 
Japan, from the end of October, 1945 to the begin- 
ning of April, 1946. It should be stressed, however, 
that over most of Asia and the Pacific the material 


is very sparse. Another trough, which was well- 
developed most of the time but showed greater 
variations in position and intensity, was situated 
generally between 90 and 50° W, i. e. near the east 
coast of North America. A third trough was indi- 
cated during much of the time at some longitude be- 
tween 0 and 90° E (Europe and western Asia); during 
the month of February, 1946, this was a major trough 
of almost the same intensity as the two mentioned 
above. 

As for the ridges, the variability appears to be some- 
what greater. The longitudes at which they occurred 
most frequently are: o—30° W (eastern Atlantic); 
100—130° W (western part of North America); 
70—90° E (western part of Central Asia). 

The major troughs and ridges (mean troughs 
as termed by Namias and Crapp (1944)) were, 
without any noticeable exeption stationary or slow- 
moving, their displacement usually being less than 
5° long. per day. Some cases of a very slow retro- 
gression of a trough or a ridge may be found. 

In several cases, individual troughs are seen mov- 
ing east within the general area of a major trough. 
The phase velocity of these individual or “daily” 
troughs, some of which are quite deep, usually 
amounts to appr. 10° long. per day. There are a 
few examples, also, of individual highs (ridges) 
moving east within the general area of ridge with a 
speed of $—10° long. per day. In'some other cases, 
troughs or ridges are seen breaking through an 
area of slight deviations, and in such cases the speed 
may amount to as much as 15° long. per day even 
for the ridges. 

A few examples of the dispersion phenomenon 
may be found in any of the months. Some of these 
examples are rather speculative, but the best ones 
(occurring in Oct. 1945 and Febr. 1946) are almost 
as striking as those in early Nov. 1945. The group 
velocities are remarkably uniform, 25—35° long. 


per day. 


5. Comparison with normal conditions. 


By the courtesy of Mr. Namias, I have had the 
opportunity to compare the results outlined above 
with a series of unpublished charts showing the mean 
height of the soo mb level (one chart for each 
month). These charts indicate, as well as those pub- 
lished by SCHERHAG (1948), that the existence of two 
main troughs, which seems to be a normal feature, 
may be interpreted as a reflection of the fact that the 
central polar low of the upper atmosphere is not of 
a circular but rather of an elliptical shape. The longi- 
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tudinal distance between the troughs is not exactly 
180° but this is in agreement with the position of 
the axis of the low: the mean maps indicate that 
this axis does not go through the north pole itself. 

The existence of a third trough is also indicated 
on some of the mean maps, the contour lines having 
a more or less triangular shape. As for the positions 
and intensity of the main troughs and ridges, there 
are some significant differences between our dia- 
grams and the normal charts. Judging from the well- 
known great anomalies of surface pressure which 


~ 


occur in middle and high latitudes, notably in win- 
ter, and from the positive correlation between 
pressure and tropospheric mean temperature which 
is known to exist for at least some areas where large 
pressure deviations are frequent (Europe and Eastern 
Atlantic), there can be no doubt that the anomalies 
of the upper air currents of individual months as 
compared with the normal conditions for that month 
are of great significance. The prediction of such 
anomalies may be called a basic problem of long- 
range weather forecasting. 
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BiO.On REVIEW 


WERNER WERENSKIOLD, Fysisk geografi I, Oslo, 
1948. 319 pp., 127 figures. H. Aschehoug & Co. 


In 1925 Mr. Werenskiold, professor of geo- 
graphy at the University of Oslo, published a 
textbook in Norwegian on Physical Geography. 
Now, after having retired from his chair, he has 
prepared a revised edition. This textbook is rather 
unique in that it treats, within the limited space 
of about 300 pages, all the main branches of geo- 
physics: size and shape of the earth, earthquakes, 
magnetic and electric state of the earth, meteoro- 
logy and oceanography. The chapter on meteoro- 
logy comprises about half of the book. 

The author has great experience in different 
kinds of geophysical field work. At the same time 
he is known as an able mathematician. His book 
contains not only purely descriptive parts, writen 
in a very concise style, but also theoretical devel- 
opments of an elementary but nevertheless exact 
character. 

The first chapter gives an outline of potential 
theory, applied to gravity, and comments upon the 
notions of isostasy and displacement of continents. 

In the second and third chapters the elements 
of seismology, terrestrial magnetism and terrestrial 
electricity are treated. In spite of the fact that Nor- 
wegian scientists have contributed fundamentally 
to the knowledge of aurora borealis, the author 
has reserved one page only for this topic. It might 
be hoped that in another edition more space might 


be given to this item, as well as to the ionosphere, 
about which so much new information has become 
available in the last few years. 

The fourth chapter gives an outline of meteoro- 
logy which seems to be well suited to make stu- 
dents of geography acquainted with the main con- 
ceptions of that branch of geophysics. The author 
presents not only the fundamental ideas concerning 
the physics of the atmosphere and climatology, but 
also the principal theoretical considerations of dy- 
namic meteorology. As to the theoretical part the 
following critical comments may be allowed. The 
derivation of the Coriolis force seems to be some- 
what more complicated than necessary. The theory 
of viscosity has been developed only for laminar 
motion. Turbulence is mentioned only incidentally 
and lateral mixing or stresses not at all. The circu- 
lation theorems and their derivatives have been 
omitted. 

The last chapter contains a descriptive and theo- 
retical treatment of oceanography (depths, proper- 
ties of the sea-water, temperature, currents, waves, 
tidal currents). The theory of the main types of 
wave motion is given in a very clear way. 

Taking into account the good qualities of this 
textbook, one feels inclined to suggest that it be 
revised and translated from Norwegian to a lan- 
guage more widely spread. So an old tradition 
might be revived: the well-known textbook of 
Mohn was once translated from Norwegian and 
made its way to many countries. 


Erik Björkdal 


Swedish Meteorological and Hydrological Institute. 
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